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Abstract 

We investigate tiie nucleon-nucleon interaction by using the meson exchange model and tlie two- 
body Dirac equations of constraint dynamics. This approach to the two-body problem has been 
successfully tested for QED and QCD relativistic bound states. An important question we wish 
to address is whether or not the two-body nucleon-nucleon scattering problem can be reasonably 
described in this approach as well. This test involves a number of related problems. First we 
must reduce our two-body Dirac equations exactly to a Schrodinger-like equation in such a way 
that allows us to use techniques to solve them already developed for Schrodinger-like systems 
in nonrelativistic quantum mechanics. Related to this we present a new derivation of Calogero's 
variable phase shift differential equation for coupled Schrodinger-like equations. Then we determine 
if the use of nine meson exchanges in our equations give a reasonable fit to the experimental 
scattering phase shifts for n — p scattering. The data involves seven angular momentum states 
including the singlet states ^5o, ^Pi, and the triplet states ^Po, ^Pi, ^Si, ^Di. Two models 
that we have tested give us a fairly good fit. The parameters obtained by fitting the n — p 
experimental scattering phase shift give a fairly good prediction for most of the p — p experimental 
scattering phase shifts examined (for the singlet states ^^o, and triplet states ^Pq, ^Pi)- Thus 
the two-body Dirac equations of constraint dynamics present us with a fit that encourages the 
exploration of a more realistic model. We outline generalizations of the meson exchange model for 
invariant potentials that may possibly improve the fit. 



* hcrater@utsi.edu 
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I. INTRODUCTION 



In this paper |Q, we obtain a semi-phenomenological relativistic potential model for 
nucleon-nucleon interactions by using two-body Dirac equations of constraint dynamics 
0, 0, H, 1^, ^ , and Yukawa's theory of meson exchange. In previous work Long and 
Crater have derived the two-body Dirac equations for all non-derivative Lorentz invari- 
ant interactions acting together or in any combinations. They also reduced the two-body 
Dirac equations to coupled Schrodinger-like equations in which the potentials appear as co- 
variant generalizations of the standard spin dependent interactions appearing in the early 
phenomenological works in this area |T^, |Tl|, |T2|, |1^, |1^, |T5| based on the nonrelativis- 



tic Schrodinger equation. This allows us to take advantage of earlier work done by other 
people on the nonrelativistic Schrodinger equation. In particular we use the variable phase 
method developed by Calogero and Degasparis JIB] , for computation of the phase shift 



from the nonrelativistic Schrodinger equation, presenting a new derivation for the case of 
coupled equations. Our potentials for different angular momentum states are constructed 
from combinations of several different meson exchanges. Furthermore, our potentials, as 
well as the whole equations, are local, yet at the same time covariant. This contrasts our 



approach with other relativistic schemes such as those by Gross and others [T^, |T^, ^ 
It is the aim of this paper to see if the meson exchanges we use are adequate to describe 
the elastic nucleon-nucleon interactions from low energy to high energy (<350 MeV) when 
using them together with two-body Dirac equations of constraint dynamics. 

Although numerous relativistic approaches have been used in the nucleon-nucleon scat- 
tering problem, none of the other approaches have been tested nonperturbatively in both 
QED and QCD as they have been with the two-body Dirac equations of constraint dynam- 
ics |, I, ^, Unlike the earher local two-body approaches of Breit |2|, |2|, 0, the 



relativistic spin corrections need not be treated only perturbatively. This means that we can 
use nonperturbative methods (numerical methods) to solve the two-body Dirac equations. 
This is a very important advantage of the constraint two-body Dirac equations (CTBDE). 
The successful numerical tests in QED and QCD gives us confidence that they may be 
appropriate relativistic equations for phase shift analysis of nucleon-nucleon scattering. 

In section II we introduce the two-body Dirac equations of constraint dynamics. In section 
III we obtain the Pauli reduction of the two-body Dirac equations to coupled Schrodinger- 



like equations. We go a step further than that achieved in the paper of Long and Crater 
in that we ehminate the first derivative terms that appear in the Schrodinger-hke equation. 
This is relatively simple for the case of uncoupled equations but not so for the case of coupled 
Schrodinger-like equations. The reason we perform this extra reduction is that the formulae 
we use for the phase shift analysis, the variable phase method developed by Calogero, has 
been worked out already for coupled equations, but ones in which the first derivative terms 
are absent. This step then becomes an important part of the formalism, allowing us take 
advantage of previous work. In section IV, we discuss the phase shift methods used in our 
numerical calculations, which include phase shift equations for uncoupled and coupled states 
and the phase shift equations with Coulomb potential. In section V we present the models 
used in our calculations, including the expressions for the scalar, vector and pseudoscalar 
interactions, and the way they enter into our two-body Dirac equations with the mesons 
used in our fits. In section VI we present the results we have acheived and in section VII 
are the summaries and conclusions of our work. 

II. REVIEW OF CONSTRAINT TWO-BODY DIRAC EQUATIONS 

The two-body Dirac equations that we will use for studying nucleon-nucleon interaction 



bear a close relation to the single particle equation proposed by Dirac in 1928 p8 



[a-p + pm + V{r)]ip = Eip. (2.1) 

For interactions that transforms as a time-component of a four vector and world scalar we 
have V{r) = A{r) + l3S{r). Of course, the single particle Dirac equation is not suitable to 
describe systems such as the mesons, (quarkonium), muonium, positronium, the deuteron 
and nucleon-nucleon scattering because the particles may have equal or near equal mass. 

The earliest attempt at putting both particles on an equal footing was in 1929 by G. Breit 
25| |2^ However, the Breit equations do not retain manifest covariant form and in QED 



the equations cannot be treated nonperturbatively beyond the Coulomb term [^, ^ . There 



have been many attempts to by pass the problems of the Breit equation and also of the full 
four dimensional Bethe-Salpeter equation. These are discussed in a number of different 
contexts in p, ^, §,0. The approach of the CTBDE provides a manifestly covariant yet 
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three dimensional detour around many of the problems that hamper the implementation 
and application of Breit's two-body Dirac equations as well as the full four dimensional 
Bethe-Salpeter equation (see also [0). In addition, as mentioned above, the approach can 
by a Pauli reduction give us a local Schrodinger-like equation. 

The CTBDE make use of Dirac's relativistic Hamiltonian formalism. In a series of papers 
(in addition to those cited above see also ||^, Crater and Van Alstine have incorpo- 
rated Todorov's effective particle idea developed in his quasipotential approach ||33| into the 



framework of Dirac's Hamiltonian constraint mechanics ||3^ for a description of two body 
systems. Their approach yields manifestly covariant coupled Dirac equations. The stan- 
dard reduction of the Breit equation to a Schrodinger-like equation for QED yields highly 
singular operators (like 6 functions and attractive potentials) that can only be treated 
perturbatively. In the treatment of the CTBDE for QED ||2^, for example, one finds 
that all the operators are quantum mechanically well defined so that one can therefore use 
nonperturbative techniques (analytic as well as numerical) to obtain solutions of bound state 
problems and scattering. (A quantum mechanically well defined potential is one no more 
singular than — l/4r^. If it is not quantum mechanically well defined, it can only be treated 
perturbatively.) Although it is encouraging that good results have been obtained for QED 
and QCD meson spectroscopy, that is no guarantee that the formalism so developed will 
lead to effective potentials in the case of nucleon-nucleon scattering that render reasonable 
fits to the phase shift data. 

Using techniques developed by Dirac to handle constraints in quantum mechanics and the 
method developed by Crater and Van Alstine, one can derive the two-body Dirac equations 
for eight nonderivative Lorentz invariant interactions acting separately or together ||35[|,P|. 
These include world scalar, four vector and pseudoscalar interactions among others. We can 
also reduce the two-body Dirac equations to coupled Schrodinger-like equations even with 
all these interactions acting together. Before we test this method in nuclear physics in the 
phase shift analysis of the nucleon-nucleon scattering problems we review highlights of the 
constraint formalism and the form of the two-body Dirac equations. 
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A. Hamiltonian Formulation Of The Two-Body Problem Prom Constraint Dy- 
namics 

Dirac extended Hamiltonian mechanics to include conjugate variables related by 
constraints of the form 0(5, p) = 0. For constraints, we may write 

(j)n{q,p)^0, n = 1,2,3, ■■,N (2.2) 
With these constraints the Hamiltonian of the system (with sum over repeated indices) 

H = qnPn - L (2.3) 
is not unique. The Dirac Hamiltonian H. includes the constraints 

n=H + \n(t)n. (2.4) 

in which H is the Legendre Hamiltonian obtained from the Lagrangian by means of a 
Legendre transformation. The may be functions of conjugate variables q's and p's. The 
equation of motion for any arbitrary function g (without explicit time dependence) of the 
conjugate variables q's and p's is then 

g=[g,H]. (2.5) 

Dirac called the conditional equality, ^ a " weak " equality meaning the constraints 0n ~ 
must not be applied before working out the Poisson brackets. Dirac called = a nonconditional 
equality or a " strong " equality. The equations of motion are 



9 = [9,'H] = [9,H + A„0„] = [g, H] + X^lg, 0„] + [g, A„]0„ ^ [g, H] + X^lg, (pn] (2.6) 
for 0„ ^ 0. 

In the two body system, we have two constraints (f)n{q,p) ~ 0, n = 1,2. For spinless 
particles they are taken to be the generalized mass shell constraints of the two particles 
0,0, namely 

ni=pl + ml + $i(x,pi,p2) ~ 0, 
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H2=p^ + m2 + $2(x,pi,P2) ~0, (2.7) 

where 

X = Xi — X2- (2.8) 

Dirac extended his idea of handhng constraints in classical mechanics to quantum me- 
chanics by replacing the classical constraints (pniliP) ~ with quantum wave equations 
'Pn{<l)P) I V') = 0) where q and p are conjugate variables. Thus the quantum forms for each 



individual particle constraint become Schrodinger-type equations |36| 

ni\i)) = Q for i = l,2. (2.9) 

The total Hamiltonian Ti from these constraints alone is 

n = XiHi + X2H2, (2.10) 

(with Xi as Lagrange multipliers). In order that each of these constraints be conserved in 
time we must have 

[n,,nm = ^^\^) = o. (2.11) 

so that 

[Hi, XiHi + A2?i2 



Xi]ni\^) + Xi[H.,n,]\^) + [H,, X2]n2\ij) + X2[n.,n2]}\^) = o. (2.12) 

Using Eq. ( |2.9| ), the above equation leads to this compatibility condition between the two 
constraints 

[7^l,7^2]|V^) = 0. (2.13) 

This condition guarantees that with the Dirac Hamiltonian, the system evolves such that 
the " motion " is constrained to the surface of the mass shell described by the constraints 
of Hi and H2 (|3^,1^ 0]). As described most recently in ||3^ this requires that 



$1 = $2 = $(xj 



(2.14) 



(a kind of relativistic Newton's third law) with the transverse coordinate defined by 



(2.15) 



and total momentum by 

P = Pi+P2. (2.16) 



To complete our review of the spinless case (|^) and establish notation we introduce the 
transverse relative momentum 

P=—Pi-—P2 (2.17) 
w w 

P-p = 0, (2.18) 
where the CM energy eigenvalue w is defined from 

{P^ + w^}\^) =0. (2.19) 

Taking the difference of the two constraints, 

{Pl-plm = -{ml-mlm, (2.20) 

we can show that the longitudinal or time-like components of the momenta in the CM system 
have the invariant forms 



Thus on these states we obtain 

[p^ + ^{xi_)-h'^{w'^,ml,ml)]\il}) = 0, (2.22) 

where 

62(w^ ml, ml) = el- ml = el -ml = {w^ - 2w^(ml + ml) + (ml - ml)H . (2.23) 

and indicates the presence of exact relativistic two-body kinematics. (By this statement we 
mean that classically — IP' = 0, would imply w = \/p'^ + ml + \/ p^ ^ ^i)- Note that 
both of the constituent invariant CM energies e\ and £2 are positive for positive total CM 
energy w greater than the square root of \m\ — . This is a direct consequence of the 
Eq.( p.20D which in turn depends on the "third law" condition necessary for compatibility . 



w + m 



mt 



2w 

w'^ + m2 — ml 



2w 



(2.21) 
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In our scattering applications below this guarantees that nucleons cannot scatter into a final 
state having an overall positive energy but with constituent positive and negative energy 
nucleons. 

In the center-of-momentum system, p = p± = (0, p), x± = (0, r) and the relative energy 
and time are removed from the problem. The equation for the relative motion is then 

{p2 + <l>(r)-62}|^) = 0, (2.24) 

which is in the form of a non-relativistic Schrodinger equation (with 2ml^ $, 2mi?jvi? ~^ 
6^). Thus the relativistic treatment of the two-body problem for spinless particles gives a 
form that has the simplicity of the ordinary non-relativistic two-body Schrodinger equation 
and yet maintains relativistic covariance. Spin and different types of interactions can be 



included in a more complete framework P, ^ ^ and will be reviewed later in this 
section. 

In addition to exact relativistic kinematical corrections, Eq. (|2.24|) displays through the 
potential $ relativistic dynamical corrections . These corrections include dependence of the 
potential on the CM energy w and on the nature of the interaction. For spinless particles 



interacting by way of a world scalar interaction S, one finds ^ 

$ = 2m^S + (2.25) 

where 

= 2.26 

w 



while for (time-like) vector interactions (described by A), one finds |^T], ^ 



$ = 2e^A - ^^ (2.27) 

where 

2 2 2 

w — mf — mi 
= ^. 2.28 

2w 

For combined space-like and time-like vector interactions (that reproduce the correct energy 
spectrum for scalar QED ||3^ ) 



$ = 2ey,A-A^ + ^VMog(l - 2A/w) + i[Vlog(l - 2A/w)f. (2.29) 

The variables m^ and (which both approach the reduced mass /i = mim2/{mi + 
m2) in the nonrelativistic limit) are called the relativistic reduced mass and energy of the 

8 



fictitious particle of relative motion. These were first introduced by Todorov in his 
quasipotential approach. Thus, in the nonrelativistic limit, $ approaches 2/1(5 + A) for 
combined interactions. In the relativistic case, the dynamical corrections to $ referred to 
above include both quadratic additions to S and A as well as CM energy dependence through 
and ew The two logarithm terms at the end of Eq.( |2.29| ) are due to the transverse or 
space-like part of the potential. Without those terms, spectral results would not agree with 
the standard (but more complex) spinless Breit and Darwin approaches (see references in 
H including ||). 

Eq.(|;24D provid es a useful way to obtain the solution of the relativistic two-body prob- 
lem for spinless particles in scalar and vector interactions and as reviewed below has been 
extended to include spin. In that case they have been found to give a very good account 
of the bound state spectrum of both light and heavy mesons using reasonable input quark 
potentials. 

These ways of putting the invariant potential functions for scalar S and vector A inter- 



actions into $ will be used in this paper for the case of two spin one-half particles (see ( |2.67 



) to ( |2.71|) and ( |4.1| ) to ( |4.3|) below). These exact forms are not unique but were motivated 
by work of Crater and Van Alstine in classical field theory and Sazdjian in quantum field 
theory Other, closely related structures will also be used. These structures play a 

crucial role in this paper since they give us a nonperturbative framework in which 5* and A 
appear in the equations we use. This structure has been successfully tested (numerically) 
in QED (positronium and muonium bound states) and is found to give excellent results 
when applied to the highly relativistic circumstances of QCD (quark model for mesons). An 
important question we wish to answer in this paper is whether such structures are also valid 
in the two-body nucleon-nucleon problem. This is an important test since the quadratic 
forms (see e.g. Eqs.( p.25D , ( p.27|) ) that appear could very well distort possible fits based on 
Yukawa-type potentials with strong couplings. 

Before going on to describe the constraints for two spin-one-half particles we mention 
an important but often overlooked aspect of the foundations of the generalized mass shell 
contraint equations given in (p.9|). It involves their derivation from an alternative starting 



point. In addition to the connection with the Bethe-Salpeter equation described in p6 



there exists a connection between constraint dynamics and Wigner's early formulation of 
relativistic quantum mechanics . In particular, Polyzou |Q has demonstrated that the 
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assumption of both Poincare invariance and manifest Lorentz covariance forces the scalar 
product for quantum mechanical state vectors to be interaction dependent. So, whereas for 
a free particle the kernal involved in the scalar product has the form 5(p^ + m^)6'(p°), in cases 
of interactions the self-adjoint nature of the kernal demands the forms 5(7ii)5(7i2) with 
compatible constraints (a related use of such delta functions to construct the state vectors 
themselves is discussed in P7|). 

B. Two spin-one-half particles 

We continue our review in this section by introducing the two-body Dirac equations of 
constraint dynamics. The Dirac equations for two free spin-one-half particles are 



5io^= (^^i-pi+mi^5i)|V^) = 

520^ = (^2 ■ P2 + m2^52) \^) = (2.30) 

where ip is the product of the two single-particle Dirac wave functions (these equations are 
equivalent to the free one-body Dirac equation). The "theta" matrices are related to the 
ordinary gamma matrices by 

0t = ^Jll5^1^, /i = 0,l,2,3, 1 = 1,2 



hi = i\j ^l5i (2.31) 
and satisfy the fundamental anticommutation relations 

[^5^,%.]+ = -l. (2.32) 

It is much more convenient to use the "theta" matrices instead of the Dirac gamma matrices 
for working out the compatibility conditions. In the reduction of complicated commutators 
to simpler form one uses reduction brackets that involve anticommutators for odd numbers 
of "theta" matrices and commutators for even numbers of "theta" matrices and coordinate 
and momentum operators. This property follows from the relation of the "theta" matrices 
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to the Grassmann variables used in the pseudoclassical form of the constraints (see 0, ^ |^). 
These fundamental anticommutation relations guarantee that the Dirac operators Sio and 
^20 are the square root of the mass shell operators —\{Pi+mD and — |(P2+"^2)- Differencing 
these implies that the relative momentum p in Eq.( p.l8D satisfies P ■ plip) = 0. 
Writing pi and p2 in terms of the total and relative momenta we obtain 



5io^ = {ei± ■ p + eiOi ■ P + rrne^iW = 

52oV' = i-02± ■ P + £2^2 ■ P + ^2^52) = 0. 



(2.33) 



The projected "theta " matrices then satisfy 



e^■p,ei■p 



J + 



0. ■ P. C 



where 



(2.34) 



(2.35) 



Defining af, = 26 1 ■ P9^,, and /3j = 29i ■ PO^i, the above two-body Dirac equations become 



(«! ■ p + f3imi)ilj = eiip 
-a2-p + I32m2)il) = €2^1) 



(2.36) 



which have the form of single free particle Dirac equations. 

Recall that in the spinless case we had the compatibility condition 



[7^1,7^2]!^) = 0. (2.37) 

It was a requirement that followed in the classical case (or the Heisenberg picture in the 
quantum case) from the individual constraints Tii being conserved in time. Similarly here 
with Si designating the form of the Dirac constraint with intereactions present, the com- 
mutator condition guaranteeing that the Dirac equations for two spin ^ particles form a 
compatible set is 

[Si,S2M = Q. (2.38) 
(It would follow from an Ti as in Eq.( |2.10D composed of a sum of the Si.) 
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We found that even for the simplest interaction, a Lorentz scalar, the naive replacement 
such as making the minimal substitutions (corresponding in the case of the single particle 
Dirac equation to Eg. with V{r) = pS{r) ), 



m, Mi(r) =mi + Si i = 1,2 (2.39) 

doe not lead to compatible constraints. Rather than detail here the earlier work steps that 
were taken to make the interactions meet the compatibility condition for scalar interac- 



tions 0, ^ we present here the form of the compatible constraints for general covariant 
interactions 

Si\^) = (cosh(A)Si + sinh(A)S2)|^/') = 

S2\^) = (cosh(A)S2 + sinh(A)Si)|^/') = (2.40) 
where the operators Si and S2 are auxiliary constraints of the form 

SilV') = (5iocosh(A) +52osinh(A))|V') = 

SslV") = (52ocosh(A) +5iosinh(A))|z^) = 0. (2.41) 
Both of these sets of constraints R][PDl PH] are compatible 



[S,,S2W = (2.42) 

[Si,S2]|^) = 0. (2.43) 

provided only that 

A(x) = A(xx). (2.44) 

Furthremore 

P-p\^) = 0, (2.45) 

the same constraint equation on the relative momentum p as in the spinless case. 

The covariant potentials are divided into two categories, four "polar" and four "axial " 
interactions. The four polar interactions (or tensors of rank 0,1,2) are 
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scalar 

Al = -1951052 = -^Ou Oi = -751752, (2.46) 

time-like vector 

Aj = JOi ■P92-P = = A/52^0i, (2.47) 

space-like vector 

Ag = ge^s. ■ 02± = = 7i± ■ 72±|0i, (2.48) 

and tensor (polar) 

= 4^0ix ■ ^2±^52^5i^i ■P92-P = = «i • a2^0i. (2.49) 

We may use each in Eqs. (|2.4CI| ) and (p.41| ) separately or as a sum 

Ap = A,, + Aj + Ag + A^ (2.50) 

to generate the sets of two-body Dirac equations with corresponding interactions. A par- 
ticularly important combination occurs for electromagnetic inteactions. While time - and 
space- like vector interactions are characterized by the respective matrices /5i/?2 and 7i_l-72±, 
a potential proportional to 71 ■ 72 would correspond to an electromagnetic-like interaction 
and would require that J = —Q. 

. iO,-0^,gi^., ^ 7La|Mo,. (2.51) 

The four "axial" interactions (or pseudotensors of rank 0,1,2) are 
pseudoscalar 

^c = ^ = Sij = -751752^(^1, (2.52) 



time-like pseudovector 



Ah = -2Hei ■ Pe2 ■ P^51^52 = -^2 Y = ^751/52752^01, (2.53) 



space-like pseudovector 



A, = -219^^ ■ 92±95i952 = -sj- = -75i7i± ■ 75272±^Oi, (2.54) 



and tensor(axial) 



Ay = -2Y9^^ ■ 92±9^ ■P92-P = -84^ = -a^ ■ (y^O^. (2.55) 
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Crater and Van Alstine found llSSi that these and their sum 



A„, = Ac + Ah + Ai + Ay (2.56) 

would be used in Eqs.( p.4(i| ) and Eqs.( |2.41"D but with the sinh(Aa) terms in Eqs.( |2.4(iD ap- 
pearing with a negative sign instead of the plus sign as is the case polar interactions . There 
is no sign change in Eqs.( f^.4l| ) for A^. 



For systems with both polar and axial interactions |^ , one uses Ap — Aa to replace A 
in Eqs. ( 12.401) , and Ap + A^ to replace the A in Eqs.(|]4l|). L, J, Q, C, E, I, Y are 
arbitrary invariant functions of x±. In this paper, we include only mesons corresponding to 
the interactions L, J, Q{J = —Q), and C. Thus we are ignoring tensor and pseudovector 
mesons, limiting ourselves to vector, scalar and pseudoscalar mesons, all having masses 
less than or about 1000 MeV. We are also ignoring possible pseudovector couplings of the 
pseudoscalar mesons. 

For computational convenience we have found it necessary to transform the Dirac equa- 
tions to "external potential" form. We obtain these forms by combining the two sets of 
equations 

Si\tp) = [cosh(A)(5iocosh(A) +52osinh(A)) + sinh(A)(52o cosh(A) + 5io sinh(A))]v:') = 



S2\ip) = [cosh(A)(52oCosh(A) +5iosinh(A)) + sinh(A)(5ioCosh(A) +52osinh(A))]?/') = 0. 

(2.57) 

and bringing the Sio operators through to the right. References [H, |3^ gives the "external 



potential " forms of the constraint two-body Dirac equations for each of the eight interaction 
matrices, A^, Aj, Ag, Ajr, Ac, Ah, A/, Ay acting alone. These forms are similar in 
appearance to individual Dirac equations for each of the particles in an external potential. 
In appeared also the form with all eight interactions acting simultaneously 

{exp{g+J^S2 + IOi + Y02)[9i-p-^e2-d{LOi-J02-gOs- m - CSi + HS2 + 1 S3 + YSi)] 
+ei cosh( JO2 + J^C4 + HS2 + Y£i)ei ■P + t2 sinh( JO2 + J^O^ + H£2 + ^^4)^2 ■ P 
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+mi cosh(-LCi + TO^ + EE2 + /^3)^5i + sinh(-LOi + TO^ + i/^a + /^3)^52}|^) = 0, 

(2.58) 

{-exp(6? + J^^2 + /Ol+m)[^2-p-^^l-9(LC'i-JC'2-6?C'3-TO-C^^l + ^^2 + /^3 + >'^4)] 

+ei sinh( JC2 + TO^ + i/^2 + YU)Oi • P + £2 cosh( JO2 + J^C4 + + ^-^4)^2 • P 



+mi sinh(-LOi + + EE<i + 1^3)^51 + ^2 cosh(-LCi + + HE^ + 1^3)^52} |^) = 0. 

(2.59) 

What is remarkable is that the above hyperbohc and exponential structures account for 
all of the "interference" terms between the various interactions. The interactions acting 
separately or in subgroupings are simple reductions of the the above. For example, in the 
case of the combined scalar, time-like, space-like and pseudoscalar interactions used in this 
paper, 

A = Aj + Ai + Ag + Ac (2.60) 
and the two-body Dirac equations (|2.58|J2.5^) reduce to 



Sx\^) (2.61) 

= (exp(e;)ei ■ V + ^1^1 ■ P + ^1^51 + I'-^^^e^ ■ d{gOs + JO2 - lOi + c^i))|^) = 



52IV) (2.62) 

= (- exp(e;)^2 ■ P + £262 ■ P + M2^52 - t^^^ei ■ d{go3 + JO2 - lOi + csi))\^) = 0. 

where 

Ml = mi cosh(L) + m2 sinh(L), 



M2 = m2 cosh(L) + mi sinh(L), (2.63) 
El = ei cosh(J) + €2 sinh(J), 
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E2 = 62 cosh( J) + ei sinh( J) . (2.64) 

In the limit mi — > oo (or m2 oo ), (when one of the particles become infinitely mas- 
sive), the extra terms dQ, dJ,dL and dC in Eqs.( p.61|p.62|) vanish, and one recovers the 
expected one-body Dirac equation in an external potential. The above two-body Dirac 
equations (without pseudoscalar interactions ) have been tested successfully in quark model 
calculations of the meson spectra 

We may rewrite the "external potential form" of the CTBDE for two relativistic spin-one- 
half particle interacting through scalar and vector potentials as (see Eqs. ( p.61| , pl62D without 
the pseudoscalar interaction) 

Silil^) = 751(71 ■ (pi - Ai) + mi + ^i)|^) = (2.65) 

S^li^) = 752(72 ■ {P2 - A2) +1712 + S2m = 0. (2.66) 

and Si introduce the interactions that the i th particle experience due to the presence 
of the other particle and are both spin-dependent order to identify these 

potentials we use Eqs.( p.61|p.62| ), and ( |2.63| , pl6^ . Then we find that the momentum de- 
pendent vector potentials A^ are given in terms of three invariant functions 010] G, -Ei, £"2 

At = ((ei - El) - ^-72 ■ -^72 • P)P + (1 - G)p^^ - -dG-^2±l^^ (2.67) 

= ((e2 - E2) - i-ii ■ -^71 ■ P)P + (1 - G)jf - -9G-7i^7rx, (2-68) 

where 

G = exp(^), (2.69) 

(with = — 1, where P = P/w) while the scalar potentials 5*^ are given in terms of three 
invariant functions [|, |[|| G, Mi,M2 

51 = Mi-mi-'-G^2-^ (2.70) 

52 = M2-m2-'-Gji-^. (2.71) 
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In QCD, the scalar potentials Si are semi-phenomenological long range interactions. The 
vector potentials are semi-phenomenological in the long range while in the short range 
are closely related to perturbative quantum field theory |Q. Of course this rewrite does 
not change the fact that Si and ^2 still satisfy the compatibility condition Eq.( p^.42D 



III. PAULI REDUCTION 



Now one can use the complete hyperbolic constraint two-body Dirac equations 
Eqs.( |2.58| , p.59| ), to derive the Schrodinger-like eigenvalue equation for the combined interac- 
tions: L{x±), J{x±), H{x±),C{x±),Q{x±), J-'{x±), I{x±),Y{x±) 0. In this paper, however, 
we include only mesons corresponding to the interactions L, J, Q{J = —Q), C, thus limiting 
ourselves to vector, scalar and pseudoscalar interactions. The basic method we use here has 
some similarities to the reduction of the single particle Dirac equation to a Schrodinger-like 
form (the Pauli-reduction) and to related work by Sazdjian [[F[P5|. 

The state vector {tp) appearing in the two-body Dirac equations ( |2.58| , p39D is a Dirac 
spinor written as 



1^) 



(3.1) 



where each is itself a four component spinor. lip) has a total of sixteen components 
and the matrices Oj's, Si^s are all sixteen by sixteen. We use the block forms of the gamma 
matrices given by Eq.(4.2) in Ref.[§ and 



1,2. 



(3.2) 



The Ef are four- vector generalizations of the Pauli matrices of particles one and two. In 
the CM frame, the time component is zero and the spatial components are the usual Pauli 
matrices for each particle. Appendix A details the procedure that leads to a second-order 
Schrodinger-like eigenvalue equation for the four component wavefunction = I'?/')! + 
in the general form 



(p' + $(r, p,o-i, 0-2, w))\<P+) = b\w)\(P+). 



(3.3) 
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Below we display all the general spin dependent structures in $(r,p,<Ti,<T2,u') explicitly, 
ones very similar to those appearing in nonrelativistic formalisms such as seen in the older 
Hamada- Johnson and Yale group models (as well as the nonrelativistic limit of Gross's 
equation). Simplification of the final result in Appendix A by using identities involving cri 
and <T2 and grouping by the term , Darwin term (r ■ p), spin-orbit angular momentum 
term L • (cri + (T2), spin-orbit angular momentum difference term L • (cti — <T2), spin-spin 
term (<Ti • <T2), tensor term (<Ti • f)(<T2 • f), additional spin dependent terms L • (<Ti x (T2) 
and (<Ti • f)(<T2 • p) + (<T2 • f)(o"i • p) and spin independent terms we obtain 

{p. _ _ ^,Mi±MiEhj + LY]r . p - lv'e-\g'^ 

L • (o-i + 0-2) _ lE2M2±MiEi _ L • (0-1 - 0-2) 1 E2M2 - MiEi 



^\ 1 2j\.^ 2 V ^ \ ^ ) 2r2 r ^ 



+ (0-1 • r)(o-2 • i)[-\v\-c + J-L)- ^v^g - g'{-c + J-L)'- g'^ + 



+^^^ p (-^ + L)' - ^^-^((o-i • r)(o-2 • p) + (^2 ■ r){cT, ■ p))}|0+) 



exp(-2e?)B'|0+). (3.4) 



where 



V = E1M2 + E2M1 



2 = Ei^ - Mi^ = E2^ - M2^ 
= h^{w) + (ei + e2)sinh^( J) + 2eie2 sinh( J) cosh( J) 
— {m\ -\- m2)sinh^(L) — 2mim2 sinh(L) cosh(L). (3.5) 
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Ei, Mi,C, J, L,Q are all functions of the invariant r. We point out that Eq.( p.4|) differs from 
the forms presented in p]. Whereas the above equation involve four component spinor wave 
functions, the ones given in |^ are obtained in terms of matrix wave functions involving one 
component scalar and three component vector wave functions. The form we choose in this 
paper is easier to compare with the earlier existing nonrelativistic forms. 

All of above equations when reduced to radial form have first derivative terms (from the 
f ■ p and (cTi ■ r)(<T2 ■ p) + (<T2 ■ • p) terms). These can be easily eliminated for the 

uncoupled equations but are problematic for the coupled equations. The variable phase 
method developed by Calogero |jl6l for computation of phases shifts starts with coupled and 
uncoupled stationary state nonrelativistic Schrodinger equations which do not include the 
first derivative terms in their radial forms. An advantage of the above for the relativistic 
case is that they are Schrodinger-like equations. Before we can apply the techniques for 
phase shift calculations which have been already developed for the Schrodinger-like system 
in nonrelativistic quantum mechanics, we must get rid of these first derivative terms. In 
terms of the above equations we seek a matrix transformation that eliminates the terms first 
order in p. 

The general form of the eigenvalue equation given in Eq.(|3.4|) is: 



g' 

[p^ - ig'r ■ p + —L ■ {ai + erg) - ih'{(Ti • f ■ p + cts ■ fcri • p) 

zr 

+ kcTi ■ <J2 + n<Ji ■ v(T2 ■ r + IL ■ {(Ti - (T2) + ijL ■ (cri x (T2) + m] 

= i3V^^|0+). (3.6) 

The m term is the spin independent part involving derivatives of the potentials. For the 
equal mass case, two terms drop out (see Eq.( p.4|) ), and the above equation becomes 

[p2 - ig'r ■ p + —L ■ (cri + erg) - ih'{(Ti ■ fcrs ■ p + cts ■ fcri ■ p) 

2r 

+ k(Ti ■ (72 + nai ■ fcr2 ■ f + m]|0+) = B^e-^^\^+). (3.7) 
We introduce the spin-dependent scale change 

10+) = exp(F + K(Ti ■ fo-2 ■ f)|^+) = {A + BcTi ■ fcrs ■ f)|^+). (3.8) 
with F, K, A, B to be determined. We find that 

p|0+) = {A + BcTi ■ r(T2 ■ f)p|^+) - tiA' + B'(Ti ■ fcra ■ f)f |^+) 

- i—[{cri - CTi ■ ff)cr2 ■ r + {(T2 - (T2 ■ ff)cri • f]|^+), (3.9) 
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and 



q' q' 

|-L ■ (o-i + 0-2)10+) = {A + BcTi ■ fo-2 ■ f)|-L ■ (0-1 + 0-2)1^+) 



g' 

+—B[2(Ti ■ 0-2 - 4iro-i ■ fo-2 ■ f f ■ p + 2ir (0-1 ■ fo-2 ■ p + 0-2 • f 0-1 ■ p) - 60-1 ■ fo-2 ■ r]\ip+) ■ (3. 10) 
We thus find that 

-tg'v ■ p|0+) = (A + BcT, ■ fo-2 • r){-ig'r ■ p)|^+) + C|^+) (3.11) 

and 

- ih'{(Ti ■ fo-2 ■ p + 0-2 ■ fo-i ■ p)|0+) (3.12) 
= {A + BcTi ■ fo-2 • r){-ih'[cri ■ fo-2 • p + 0-2 ■ fo-i ■ p])|'0+) + Dlip^) 

and finally 

p'|0+) = {A + B(Ti ■ fo-2 ■ f)p'|^+) - 2i{A' + E'o-i • fo-2 ■ f)f ■ p|^+) (3.13) 
2B 

+ i— [2(Ti ■ fo-2 • f f ■ p - (0-1 ■ fo-2 ■ p + 0-2 ■ f 0-1 • p] 1^/'+) + E\iIj^) 

where C and D and E do not involve p and are given by 

C = -(7'(A' + 5'o-i-fo-2-f), (3.14) 



T3 

D = -2h'i(Ti ■ fo-2 ■ f A' + B') - 2h' — \L ■ (0-1 + 0-2) + 2 - 0-1 • fo-2 ■ f + 0-1 • 0-2], (3.15) 

r 

and 

2 B 

E = -{A" + B"ai ■ fo-2 ■ f) - -(A' + B'ai ■ fo-2 • f) - 2— (0-1 • 0-2 - 3o-i ■ fo-2 ■ f). (3.16) 

The general form of the eigenvalue equation then becomes after some detail |jl| 

[A + Bai ■ fo-2 ■ f ) [p^ - ig'i ■ p + ■ (0-1 + 0-2) - ih'{(Ti ■ fo-2 ■ p + 0-2 ■ f 0-1 ■ p)] |V^+) 

+ i^B\2(Ti ■ 0-2 - 4zro-i ■ fo-2 ■ f f ■ p + 2zr(o-i • fo-2 • p + 0-2 ■ fo-i • p) - 60-1 ■ fo-2 ■ f] 

2r 

2B 

- 2i{A' + B'cTi ■ fo-2 ■ f )f ■ p + i— [2o-i ■ fo-2 ■ f f ■ p - (0-1 ■ fo-2 ■ p + 0-2 ■ f 0-1 ■ p)] 
+ {kcTi ■ 0-2 + ncTi ■ fo-2 ■ r){A + Bcri ■ fo-2 ■ f ) + R + m)\ip+) 

= exp(-2^)(v4 + B(Ti ■ fo-2 • f)|^/'+) (3.17) 
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in which R = C + D + E. 

Now, to bring this equation to the desired Schrodinger-hke form with no hnear p term 
we multiply both sides by 

^ ^ -.N 1 (A — B(Ti ■ rcTo ■ r) , , 

{A + BcT, ■ va, ■ r)-' = ^ ^2 _ ^^'^^^ 

and find, using the exponential form above that appears in Eq.( p.8|) , (and some detail 

{A + B(Ti ■ ■ r)-^[-2i{A' + E'c^ ■ icr^ ■ r)]r ■ p 
= -2i{F' + K'(Ti-r(T2-r)r-p, (3.19) 

and 

{A + BcTi ■ i(T2 ■ Yy^i—\lai ■ fcrg ■ ff ■ p - (cti ■ r(T2 ■ p + (T2 ■ fcri ■ p)] 

2i sinh(/r) cosh(i^') ^ / - ^ 

= [2(Ti ■ Y(T2 • rr • p - (CTi ■ VCT2 ■ p + cr2 • rcTi ■ p)J 

+ G (3.20) 
where ([0]) 

G=-^^i!^^L. + (3.21) 



and 



(A + BcTi ■ T(T2 ■ f) ^—B\2(Ti ■ (T2 — 4:irai ■ ra2 ■ f f ■ p 

2r 

+ 2ir{cri ■ y(T2 • p + 0-2 ■ i<Ti ■ p) - Q<Ji ■ fcra ■ f] 

ig' smh.(K) coshiK) , ^ , ^ ^ 

= -4rcri ■ rcTa • r r ■ p + 2r{cri ■ rcrs ■ p + 0-2 ■ rcri ■ p) 

2r 

- 2i(Ti ■ (T2 + Qi(Ti ■YCT2-r\ + H (3.22) 



where ([0) 



fif' sinh^(if) .^^ , X „ ^ „ .^ 

H = 2L ■ (o-i + (T2) - 2cri • f o-g ■ f + 2(Ti ■ CT2 + 4 . 

2r 



Note that G and H do not contain linear p type of terms. Now collect the three different 
linear p type of terms in equation (|3.17|) : 



(_2zF' - zg')r ■ p, (3.23) 
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,_^. smh{K) cosh(ir) _ ^ sinh(ir) cosh(ir))(cri ■ f 0-2 ■ p + 0-2 ■ fo"! ■ p), (3.24) 



Ui ^'^'^HK) cosh(ir) _ 2^ gj^j^.^^ cosh(ir)(7' - 2iK')(Ti ■ fcr2 ■ f f ■ p. (3.25) 

r 

If we set the first of the above equations to 0, we obtain the expected result (for the uncoupled 
portion of the equation) 

F' = -g'/2. (3.26) 

If we set h' = —K' and use p = f(f.p) — then the two expressions ( P . 24|) and (p.25|) 
combine to 



, sinh(i^) cosh(i^) , , , , , , cti ■ fcr2 ■ fL ■ (cti + cto) 

(2 ^ + /i'-^'sinh(ir)cosh(Js:))^ = r 

which contains no f ■ p. Thus the matrix scale change 

10+) = exp(-(7/2) exp(- ■ ra^ ■ f ) |^+) (3.28) 

eliminates the linear p terms. 
Further note that 

[A + BcTi ■ i(j2 ■ r)^^{kcri ■ cr^ + ncr^ ■ r(T2 ■ r){A + B<Ji ■ y(T2 ■ f) 
= {kai ■ <J2 + na-i ■ r(T2 ■ f), 

{A + BcT, ■ v(T2 ■ f)-^C|^+) = -g\F' + K'a, ■ v(T2 ■ f)|V+), (3.29) 
and (after some algebraic detail |1|]) 

{A + BcTi ■ ra-2 ■ r)-'D\ij+) = -2h\K' + F'cn ■ v(T2 ■ f)|V^+) 



, ,cosh(i^') sinh(i^') _ , n „ ^ ^ n, , v 

-2h' — ^[L ■ (<Ti + 0-2) + 2 - CTi ■ fo-2 • f + CTi ■ 0-2]|^+) 



,sinh^(i^) 

+2h [cTi ■ Ycr2 ■ rL • [cti + (T2) + 3cri ■ rcrg ■ r - cti ■ cr2\\yj+). (3.30) 
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also 



[A + Bai ■ fo-2 ■ r)-^E\ilj+) = -[F" + F'^ + K'^ + {2F'K' + K")ai ■ fcra ■ f] 



2r^/ ^ coslifi^') sinh(_ft') , 

sinli^(JC) , , , , 

+ 2 (CTI ■ fo-2 ■ f - CTi ■ 0-2 - 2). (3.31) 

So combining all terms and grouping by term , spin independent terms, spin-orbit angular 
momentum term L- (cri + cr2), spin-spin term (cr^-cr2), tensor term (<Ti-r)(cr2 ■ r), additional 
spin independent term we have our Sclirodinger-like equation 



r 

2 clr,V,2 

_ F" - F'^ - K'^ - -F' - A 



^^g\ g'sinh^K) 2sinh^(ir) ^^, cosh(ir) sinh(ir) ^ 
+ L ■ (o-i + cr2)[— + 2h J 

^ , ,,sinh2m sinh(K) coshm h' g' smh(K) cosh(K) , 
+ o-i ■ fo-2 ■ fL ■ o-i + cr2 2/i' ^ + 2 ^ + - - ^ 

(?' cosh(ii') sinh(ii') (7' sinh^(K) ,cosh(ii') sinh(i^) 
+ (Ti ■ cr2[k -\ 1 2h 

/y /yt /y 

,smh^{K) cosh(is:) sinh(ir) ?,m\i^{K) 

2/i 2 T 2 r |-|- 

3(7' cosh(i^) sinh(ii') g' smh^ K , . . . 2/i' cosh sinh 
(T\ ■ r(T2 ■ r n g K — 2h F H 

ry ry ry 

+ 6ft'!l^ - (2F'K' + A-) - Ha" + 6^°^h(/0 Binh(A-) ^ ^stah^ ^ ,^ 

^ ^ .y^ ^ ry ^ 

= B^e-^^\tlj+) (3.32) 
Comparing Eg. ( |3.7] ) with Eg. ( |3.4| ) we find 



g' = 2g' - E^Ml±MlEl(^j + L)' = 2g' - log' V = -2F', (3.33) 



h' = i^— ^ = -K', (3.34) 
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n = + J-L)- ^v'g - Q\-c + j-Ly- + ^g' 



+ + J-Ly + l \og'V{g -C + J- L)', (3.36) 

2r 2 



m = -Iv^g-lg'^ _ 1(C7 + J - L)'(-C + J - D' + ^g'log'V. (3.37) 



Eq.( p.32[) and it's derivation is an important part of this paper. It will provide us with 



a way to derive phase shift equations using work by other authors who developed methods 
for the nonrelativistic Schrodinger equation. First we need the radial form of the coordinate 
space form of this equation. 

The following are the radial eigenvalue equations for singlet states ^5*0, ^Pi, and 



triplet states Pq, -Pi, Si, Di corresponding to Eq. (|3.32|) with the above substitutions. 



We emphasis that unlike the potentials used by Reid, Hamada- Johnson and the Yale group 



jl^ , 0, |T5[, our potentials are fixed by the structures of the relativistic two-body Dirac 
equations and we do not have the freedom of choosing different potentials for different 
angular momentum states. 

^Sq, ^Pi, ( a general singlet ^Jj): For these states L ■ (cti + 0-2) = 0, cti ■ 0-2 = —3, 
cTi • ia2 ■ f = — 1. There is no off diagonal term. We find (adding and subtracting the 6^ 
term) 



where our effective potential for above equation is 

^ , , (26? - log(P) - J + L)'2 {2g - log(P) - J + L)" , {2g - log(D) - J + L)' 
•^(n = A + ^ + Z 



+^V\-C+J-L-3g)-^{C+J-L-g+2lnD)\-C+J-L-3g)'-B^e-^^+b\w). (3.38) 

Our radial eigenvalue equations for singlet states ^Sq, ^Pi, ^D2 have the same potential 
forms except for the ^^^^^ angular momentum barrier term. Later, we shall show that their 
potentials actually are different due to the inclusion of isospin ti-T2 terms. 

^Pi{ a general triplet ^Jj)'- For these, L ■ (cr^ + a2) = —2, (Ti ■ (T2 = I, cri ■ i(T2 ■ f =1 

For the ^Pi state the radial eigenvalue equation is 
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with 



(2g - log(p) + J - Ly^ , (2g-\og(v) + j-L)" , (g + j-L-cy 

5>(r) = \ \ 



-^V''{-C+J-L+g)+^{2\og{V)-{C+J-L+3g))\j-L-C+gy-B''eM-'^S)+b'{w), 

(3.39) 

The ^5*1 and "^Di are coupled states described by u^{r) and M_|_(r) and their radial eigen- 
value equations are (using L ■ (cti + (T2) = 2(j — 1), cti • (T2 = 1, <Ti • icr2 ■ f =2jTi (diagonal 
term), and (Ti • fcr2 • f =^^^^^i^(off diagonal term)) in the form 

(f 

{-T^ + $ii(r)K + ^i2{r)u+ = h'^u- (3.40) 
{-T^+^ + *22(r)}ii+ + $2i(r)ii- = 6V (3.41) 

where 

8(2^'-log'P)sinh2(/i) 8(J-L)'cosh(/i)sinh(/i) lQsmh^{h) {2g' -\o^{V)f 
^ii{r) — {- h ^ \ ^ 

, {J-Ly' , (2V-log-(P))(J-L)- (2g---log"(P)) 
4 6 2 

, {J-Lf ^ (2g--l()g-(P)) ^ (J-L)- 
6 r 3r 

i[_lv2(-c + J - L + ^?) - ^?'( J - L - c + ^?)' + ^ iog'(p)(e? + J - L - cyy 

]g'^ - \g'C - 1(C + J - L)'(-C + J - L)' - -B^ exp(-2^;) + 62(w)} (3.42) 

^ 



$^2(r) = ^{(2^' - log^g)( ^^"^^^^^ sinh(/i) _ sinh^^ + ( J - L)'(^^^^^i^^ - ^"^^(/i) sinh(/i) . 
_ 6cosh(/i) sinh(/i) ^ 2sinh^(/i) _ 1^2/^ + j _ ^ + ^) 

j,2 ^2 2 

_g^(j_L-C + g)^+ ^^^ + '^~^~^^V ^log^(P)(g + J-L-C)^ 

^ (2g--log-(P))(J-L)- ^ (J-V ^ (J-L)\ ^3 43^ 
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8 {2g' - log V) sinh2(/i) 8 (J - L)' cosh(/i) sinh(/i) 16 sinii^ih) {20' - log{V)y 



*22(r) = {--^^^ ^-^ ^ - '—-^ ^ + + 



3 7* 3 7* 3 4 
(J-Lf (2j^-log-(P))(J-Ly (2g--''- log--(P)) (J-Lr 

4 6 2 6 

2(2g^ - log^(p)) ^ 2{j-Ly {g + j-L-cy 

r 3r r 

- li'l^^-c + J-L + g)- g'{j -L-c + gy + ^ iog'{v){g + cyy 

]g'^ - ^g'C -]{C + J- Ly{-C + J-Ly-B^ exp(-2e?) + b\w)} (3.44) 



$21 (r) = $12 (3.45) 

4:\/2\(J Ll'f^^"^^^^^ + ^) 2 cosh(/i) sinh(/i) ^ (2^' - log' P) cosh(/i) sinh(/i) 
r 2r r"^ r 

(Note that because of the spin-orbit-tensor term, the potential is not symmetric.) In 
Appendix B we give the coupled equations for triplet ^jj-i and ^jj+i for general j. The 
remaining special case is for the ^Pq state and has the form 



f_ 2 
dr 

where 



^ (2g - log(P) - J + Ly^ ^ (2g - log(P) - J + Lf ^ (log(P) - (4g + J - L - 2C))' 
4 2 r 

+^V2(-C + J - L + 6?) - + ^(c" - (J - L)'2) + g'{^g ^j-L-cy 

-\ log'(P)(J - L - C + e;)' - -B" exp(-2^?) + \?{w), (3.46) 
Now we can apply the techniques already developed for the radial Schrodinger equation 

(-^ + ^^^^ + 2mV^,,(r))^ ^ 2mEv (3.47) 
in nonrelativistic quantum mechanics to the above radial equations by the substitutions 
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By comparing $ and 2m one could determine whether our $ is similar to standard type 
of phenomenological potentials such as Reid's potentials. But first, in the next section, we 
discuss the models we used in our calculation. This includes how we choose the Q, L and C 
invariant potential functions, the mesons we used in our calculation, and the way they enter 
into the two-body Dirac equations. 

IV. THE INVARIANT INTERACTION FUNCTIONS 
A. The Q and L Interaction Functions 

Our dynamics depends on how we parametrize the invariant interaction functions Q, L and 
C. We first consider how to model Q and L, corresponding to vector and scalar interactions. 
As we have seen, in order that Eg. ( p.65| ) and Eq.( p.66D satisfy Eq.( p.42| ), it is necessary that 



the invariant functions G, Ei,E2, Mi and M2 depend on the relative separation, x = X1—X2, 
only through the space-like coordinate four vector = + P^{P ■ x), perpendicular to 
the total four-momentum P. For QCD and QED applications, G, Ei,E2 are functions 010] 
of an invariant A. The explicit forms for functions Ei,E2, G are 

E, = G{e^-A) (4.1) 
E2 = G{e2 - A) 

and 

= jY^y (4-2) 

The function A{r) is responsible for the covariant electromagnetic-like A'^. Even though 
the dependencies of Ei, E2,G on A is not unique, they are constrained by the requirement 
that they yield an effective Hamiltonian with the correct nonrelativistic and semi-relativistic 



limits (classical and quantum mechanical ||40| , |4l[|). For QCD and QED application , Mi 



and M2 are functions of two invariant functions A{r) and S{r) 

M^{A, S)=ml + G\2m^S + S^) 



Ml{A,S) = ml + G\2m^S + S^). (4.3) 
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The invariant function S{r) is responsible for the scalar potential since Si = 0, if S{r) = 0, 
while A{r) contributes to the Si (if S{r) 7^ ) as well as to the vector potential A^. So, 
finally, the five invariant functions G, Ei,E2, Mi and M2 (or Q = —J,L) depend on two 
independent invariant potential functions S and A. (Compare also the spin independent 
portions to Eqs.(g]^g]^ through calculation of Ef - Mf - b^.) 

Expressing G, Ei,E2, Mi and M2 in terms of S and A is important for semi- 
phenomenological and other applications that emphasize the relationship of the interactions 
to effective external potentials of the two associated one-body problems. However, the five 
invariants G, Ei,E2, Mi and M2 can also be expressed in the hyperbolic representation 
|35| in terms of the three invariants L, J and Q (see Eqs.( p.63 ), ( 2.64 ) and ( p.69 )). L, J 
and Q generate scalar, time-like vector and space-like vector interactions respectively and 
enter into our Dirac equations via the sum + Aj + Ag where Eqs( p.46' , |2.47 , 2.48 ) define 

AL,Aj,Ag. 



We may use Eqs.( 2.41 ) to relate the matrix potentials A to a given field theoretical 
or semi-phenomenological Feynman amplitude. As mentioned earlier, a matrix amplitude 
proportional to 7^ ■ 72^ corresponding to an electromagnetic-like interaction would require 
J = —Q. Matrix amplitude proportional to either /1/2 or 71 ■P72--P would correspond to 
semi-phenomenological scalar or time-like vector interactions. The two-body Dirac equations 
in the hyperbolic form of Eq.( |2.4lD give a simple version for the norm of the sixteen 
component Dirac spinor. The two-body Dirac equations in " external potential " form, 
Eq.( p.65|) and Eq. (p.66|) , (or more generally ( p.61| ) and ( p.62| )) are simpler to reduce to 
the Schrodinger-like form and are useful for numerical calculations (see Sazdjian |^ for a 
related reduction). We describe the parametrization of the pseudoscalar interaction C below 
in Eq.(O). 



B. Mesons Used in the Phase Shift Calculations 



We obtain our semi-phenomenological potentials for two nucleon interactions by incor- 
porating the meson exchange model and the two-body Dirac equations. Because the pion 
is the lightest meson, its exchange is associated with the longest range nuclear force. The 
shortest range behaviors of our semi-phenomenological potentials are modified by the form 
factors, which are treated purely phenomenologically. We exclude heavy mesons that medi- 
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ate the ranges shorter than that modified by the form factors. The intermediate range part 
of our semi-phenomenological potentials comes from exchange of mesons which are heavier 
than the pion. We use a total of 9 mesons in our fits. These include scalar mesons cr, 
and /o, vector mesons p, uj and 0, and pseudoscalar mesons vr, r] and rj' . In this paper, we 
are ignoring tensor and pseudovector interactions, limiting ourselves to vector, scalar and 
pseudoscalar interactions, all with masses less than about 1000 MeV. See the Table below 
for the detailed features of the mesons we used |^ . 



TABLE I: Data On Mesons(T=isospin, G=G-parity, J=spin, vr =parity) 
Particles Mass(MeV) .r Width(MeV) 



7r± 


139.57018±0.00035 1" 


0- 






134.9766 ± 0.0006 


1- 


0- 






547.3 ± 0.12 


0+ 


0- 


(1.18 ±0.11) X 10 


P 


769.3 ± 0.8 


1+ 


1- 


150.2 ± 0.8 


UJ 


782.57 ± 0.12 


0- 


1- 


8.44 ± 0.09 


V 


957.78 ± 0.14 


0+ 


0- 


0.202 ± 0.016 




1019.417± 0.014 


0- 


1- 


4.458 ± 0.032 


/o 


980 ± 10 


0+ 


0+ 


40 to 100 


do 


984.8 ± 1.4 


1- 


0+ 


50 to 100 


a 


500-700 


0+ 


0+ 


600 to 1000 



C. Modeling the Invariant Interaction Functions 



We initially assume the following introduction of scalar interactions into two-body Dirac 
equations (see Eqs.(|2T0|, |2Tl|), ^) ): 

S = -sf— - (rvr.)£:— g}:-^ (4.4) 

where g^, gl^, gj^ are coupling constants for the a, Oq and /o mesons and m^, niaQ and mj^ 
the corresponding masses. {ti-T2) is 1 or —3 for isospin triplet or singlet states. 

Pseudoscalar interactions are assumed to enter into two-body Dirac equations in the form 
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(see Eq.(pl|)) 

C = {rrr.f-^- + + , (4.5) 

w r w r w r 

where w = ei + 62 is total energy of two nucleon system, gl, g'^, g"^, are coupling constants 
for mesons vr, 77 and 77' respectively and m,r, ""^r? and m^/ the corresponding masses. This 
form for C yields the correct limit at low energy. 

We also initially assume that our vector interactions enter into two-body Dirac equations 
in the form (see Eqs. ( [4.1| ) and ( |4.2| )) 

p-mpT p-irijnr p-m^r 

where (yf^, (7^, g"^ are coupling constants for mesons p, u and (p and rrip, rrio; and m</, are the 
corresponding masses. 

We use form factors to modify the small r behaviors in C and A, that is, the shortest 
range part of nucleon-nucleon interaction. We choose our form factors by replacing r in S, 
C and A with 



r2 + tI. (4.7) 

In our first model, we just use two different r'^s to fit the experimental data, one tq for the 
pion, one for all the other 8 mesons whose masses are heavier than pion's mass. We set 
these two r'^s as two free parameters in our fit. These form factors are different from the 
conventional choices, usually given in momentum space, but the effects are similar. 

In the constraint equations, A and S are relativistic invariant functions of the invariant 
separation r = \/x\_ (see below for the distinction between A and A). Since it is possible that 
A and S as identified from the nonrelativistic limit can take on large positive and negative 
values, it is necessary to modify G, Ei,E2, Mi and M2 so that the interaction functions 



remain real when A become large and repulsive These modifications are not unique 

but must maintain correct limits. 

We have tested several models, two of which can give us fair to good fit to the experimental 
data. 

a. Model 1 For Ei = G{ei — A) to be real, we need only require that G be real or ^ < 



w/2. This restriction on A is enough to ensure that Mj = G^Jmf{l — 2A/w) + 2myjS + S"^ 
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be real as well (as long as S* > ). In order that A be so restricted we choose to redefine it 

as 



A = A, 



A < 



(4.8) 



A 



A 



Vaa^T^' 



A>0. 



(4.9) 



This parametrization gives an A that is continuous through its second derivative. 

We next consider the problems that may arise in the limit when one of the masses becomes 
very large . Even though both our masses used in this paper are equal, we demand that 
our equations display correct limits. We must modify Mi and M2 so that it has the correct 
static limit (say m2 00). It does appear that Mi mi + S when 777,2 —>■ 00. However 
this is only true if mi + S* > 0. In other words, in the limit m2 — > 00, the two-body Dirac 
equations would reduce to 

(7-j9i + |mi + 5|)|V^) = 0. (4.10) 

This would deviate from the standard one-body Dirac equation in the region of strong 
attractive scalar potential {S < —mi). In order to correct this problem, we take advantage 
of the hyperbolic parametrization. We desire a form for Mj that has the expected behavior 
{Mi rrii + S in the limit when S becomes large and negative and one of the masses is 



large). So we modify our L in the following way |24 



smh L = 1 1 + 



A)S. 



w 



and for 



5 > 

Ml = ml + G^{2mu,S + S^) 



5 < 



(4.11) 



Ml = ml + G^{2m^S + S^) (4.12) 

with Eqs.( ^) . 

A crucial feature of this sinhL extrapolation is that for fixed S", the static limit(m2 ^ mi) 
form is sinhL —>■ S/w which leads to Mi mi + S. The above modifications are not unique 



give the correct semirelativistic limits 
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b. Model 2 This model comes from the work of H. Sazdjian Using a special tech- 
niques of amplitude summation, he is able to sum an infinite number of Feynman diagrams 
(of the ladder and cross ladder variety). For the vector interactions, he obtained results that 
correspond to Eq.{^^ to Eq.{^^ and Eq.(|]T]) to Eq.(^ (modified here in Eq.(|;^ for 
A > 0). For scalar interactions {L{S,A)) he obtained two results. One again agrees with 
Eq.( p2.25D and Eqs.( [4.3D . As we have seen above this must be modified (see Eq.( [4.1lD ) for 
S < 0. His second result is the one we use here for our second model for {L{S,A)). That 
replaces Eq.( 4.11 )and Eq.( 4.12|) with the model: 



then 



while if 



we let 



In both case we let 



S + A>0 



S -A+^M±£^ (4.13) 

^A{S + Ay + w^ 



S + A<0 



S — >-A + S + A. (4.14) 



sinhL = sinh(--ln(l- ^^'^^'^^ -g). (4.15) 
2 w 



D. Non-minimal Coupling Of Vector Mesons 

The coupling of the vector mesons in Eq.( f4.6| ) corresponds in quantum field theory to 
the minimal coupling gpV^ip'y'^ ip analogous to eAf^ip'y^ ip in QED. In our model, we are not 
concerned about renormalization, since the quantum field theory is not fundamental, so that 
we cannot rule out the non-minimal coupling of the p, u, analogous to 

^^i^b'.ili^F,,- (4.16) 

We can convert the above expressions to something simpler by integration by parts and 
using the free Dirac equation for the spinor field. This nonrenormalizable interaction be- 
comes 

p ApTfi AT 2p 

^i^^b'^ l^W^^ ^ -^^i^lyA, - i-{^dy - {d^i^)i^)A,. (4.17) 
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The first term can be absorbed into the standard minimal coupling while the second term 
gives rise to an amplitude written below. Changing from photon to vector mesons (p) and 
using on shell features we find 

^f,{v,. + '-^){p + PTiP + pT 4pJp + v'f -m^rn% + q^) 



M2 (g2 + m2 - le) M2 (g2 + ^2 - ie) {q^ + ml - ie) 



(4.18) 



where q = p—p'. The mass M is a mass scale for the interaction, rriN is the fermion (nucleon) 
mass and rup is the p meson mass. 

How does this interaction modify our Dirac equations? Which of the 8 or so invariants 
are affected (see Eq. (|2.46| ) to Eq.( p.56D )? In terms of its matrix structure, the above would 
appear to contribute to what we called (see Eq.( p.46| )). It is as if we include an additional 



scalar interaction with an exchanged mass of a p and subtract from it the Laplacian (the 
terms in Eq.( ^.18D ). That is 



S -^S + S' - V^S'/Aml (4. 19) 

where 

^ " M2 r ^^-^^^ 

so that the modification is rather simple. It has the opposite sign as the vector interaction. 
That is, it would produce an attractive interaction for pp scattering. But to lowest order, its 
attractive effects are canceled by the contribution of the first term on the right hand side of 
Eq.([4.17D. In our application, this means that Eq.(|4.4]) and Eq.(|4.tj|) are replaced (including 



the ro by Eq.(|3)) by 

S = -9a^- - 9ao^ 9fo— S (4.21) 

where 

S' ^ (, . - + ,2(1 - + - (4.22) 



and 

p-mpf p-ruwr p-m^r 

A = in- r,){gl + g';f-^ + {gl + g'^f-^ + (^ + g'^ f-^ (4.23) 
where g'^, g'^, g'^ are also couphng constants we will fit. 
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V. VARIABLE PHASE APPROACH FOR CALCULATING PHASE SHIFTS 



In this section, we discuss and review the phase shift methods which we used in our nu- 
merical calculations, which include phase shift equations for uncoupled and coupled states 
and the phase shift equations with Coulomb potentials. The variable phase approach devel- 
oped by Calogero has several advantages over the traditional approach. In the traditional 
approach, one integrates the radial Schrodinger equation from the origin to the asymptotic 
region where the potential is negligible, and then compares the phase of the radial wave 
function with that of a free wave and thus obtain the phase shift. In the variable phase 
approach we need only integrate a first order non-linear differential equation from the origin 
to the asymptotic region, thereby obtaining directly the value of the scattering phase shift. 

This method is very convenient for us since we can reduce our two-body Dirac equations 
to a Schrodinger-like form for which the variable phase approach was developed. Thus, 
we can conveniently use this variable phase method to compute the phase shift for our 
relativistic two-body equations. 



A. Phase Shift Equation For Uncoupled Schrodinger Equation 

Ref. [0 gives a derivation of a nonlinear equation for the phase shift for the scattering 



on a spherically symmetrical potential with the boundary condition 



ui{0) = (5.1) 
of the radial uncoupled Schrodinger equation 

<(0 + [k' - ^^^^ - V{r)]ui{r) = 0. (5.2) 

The radial wave function is real, and it defines the "scattering phase shift" Si through the 
comparison of its asymptotic behavior with that of the sine function: 

ui{r) '^-^^ const ■ sin(/cr — tt + ^i) (5.3) 
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The equation that Calogero derives is 



t'lir) = -^V{r)[ji{kr) -tiir)hi{kr)Y 



(5.4) 



where ti{r) has the hmiting value tan^^ with the boundary condition ti{0) = 0. This is 



a first-order nonhnear differential equation and can be rewritten [|16[ in terms of another 
function 6i{r) defined by 

ti{r)= tan 5i{r) (5.5) 



with the boundary condition 



and limiting value 



Si{r) ^ 



(5.6) 



lim 5i(r) = 5,(00) = 5, 



The differential equation for 6i{r) is ||T6| 



(5.7) 



6'i{r) = —k ^V{r) [cos 6i{r)ji{kr) — sin 6i{r)ni{kr)y 



(5.8) 



The solution of this first order nonlinear differential equation yields asymptotically the value 
of the scattering phase shift. The function 6i{r) is named the "phase function " and Eq. ( ^.8|) 
is called the "phase equation " . It is our main tool for studying the properties of scattering 
phase shifts. Eq.(|5.8|) becomes particular simple in the case of S waves 



S'^(r) = -k-W{r) sin^[kr + 5o{r)]. 
Now, since our Schrodinger-like equations in CM system has the form 



(5.9) 



[V2-62_$]|^) = (5.10) 

we can directly follow the above steps to obtain the phase shift by swapping k ^ b , and 
V ^ ^. There is no change in the phase shift equation, even though our quasipotential $ 
depends on the CM system energy w. 

We have found in convenient to put all the angular momentum barrier terms in the 
potentials, and change all the phase shift equations to the form of S state-like phase shift 
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equations |T6[ . This puts our phase shift equations in a much simpler form. For spin singlet 
states, our phase shift equations become just 



S[{r) = -r^$i(r) sin^lbr + 6i{r)]. (5.11) 

This equation is similar to the ^5*0 state phase equation(see Eq.( ^.9|) ), but it works well for 
all the singlet states when the angular momentum barrier term (^^^^) is included in $/(r), 

$,(r) = $(0 + (5.12) 

Because the nucleon-nucleon interactions are short range, we integrate our phase shift 
equations (for both the singlet and triplet states) to a distance (for example 6 fermis) where 
the nucleon-nucleon potential becomes very weak. Then the angular momentum barrier 
terms ''^^^2^^ dominate the potential ^i{r) and we let our potential $/(r) = ^^^^^^ and integrate 
our phase shift equations from 6 fm to infinity to get our phase shift. (This can be done 
analytically in the case of the uncoupled equations |]I6[ .) 



Because of the modification of our phase shift equations, we also need to modify our 
boundary conditions for phase shift equations. For the uncoupled singlet states ^Pi, 
and triplet states ^Pq, ^Pi, the modified boundary conditions are [l^ 



m = -j^b. (5.13) 

This is implemented numerically by an additional boundary conditions at r = h, so our 
boundary conditions for uncoupled singlet states ^Pi, and triplet states ^Pq, ^Pi are 

Si{h) = -j^bh (5.14) 



where h is stepsize in our calculation, b = vfo^ and of course (5/(0) = 0. So for P and D 
states, the new boundary conditions are 6i{h) = —^bh, 62(h) = —\bh respectively. 

B. Phase Shift Equation For Coupled Schrodinger Equations 

For coupled Schrodinger-like equations, the phase shift equation involves coupled phase 
shift functions. We discuss an approach here different from that originally presented in 



T^. The key idea for this new derivation is taken from one presented in a well known 
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quantum text |^ . We present an appropriate adaptation of this idea here in the uncoupled 
case to demonstrate the general idea and then extend it to the coupled case. Consider a 
radial equation of the form 

(P 

+ ^i{r))u = h^u 



dr 



Following 1^ we assume 



u{r) = A{r) sin(6r + 6i{r)) (5.15) 

■u'(r) = bA{r) cos(6r + 5;(r)). (5.16) 

Taking the derivative of the first equation we find that 

A' = -A6'iCot{br + 6) (5.17) 

and then using this and Eq.( |5.1(j| ) the above radial Schrodinger equation reduces to Eq. ( p.ll| ). 
The coupled radial Schrodinger equation has the form 

U" = -b^U + ^{<^lU + U<^l) (5.18) 

where both U and $l are two by two matrices. The effective quasipotential matrix is of 
the form 

$L = " . \ (5.19) 
while the matrix wave function is assumed to be of the form 

U = ^[^sin(6r + I?)+sin(6r + 1))^] (5.20) 

U' = ^[^cos(6r + V) + cos(6r + V)A] (5.21) 

with (using Pauli matrices to designate the matrix structure) 

V=5 + D ■ <T 

A=a + A • cr. (5.22) 

(The functions V, A are not related to earlier functions which use the same symbols). We 
further assume (for real and symmetric potentials) that both the phase and amplituded 
functions are diagonalized by the same orthogonal matrix 

U=RUR^^ = {a + Aag) sin(6r + 5 + Da^). (5.23) 
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Combining Eqs. (|5.2(3|J5.21| ) together with Eg. ( p. 181) so as to produce the analogue of the 



phase shift Eq.( |5.17D requires we use the following properties of the orthogonal matrix R 

R = 



cose[r) sme[r) 
-sme(r) cose(r) 



R'R-^ = £' \ ° ^ I = e'i(72 (5.24) 

In appendix C we derive the coupled phase shift equations {6 = {6i + 62)/2, D = {6i — 62)/2) 
below: 



S[{r) = -^[($11 + ^i^^lf^)cos2e(r) (5.25) 




+ (<f22 + ^^^^^)sin2e(r) + $ 12 sin 25 (r)] sin [hr + 5i{r)) 



^2(0 = -\[{^22 + Mfi±l)) cos2e(r) (5.26) 




($11 + M^lJlll) sin2 s{r) - $12 sin 2e{r)\ sin^(6r + 62{r)) 



1 rlr hik + 1) ^2(^2 + l)n 

^ W = r-TTT^-FTTTi^ "^11 + 2 - ^^^ 2 sin2e r 
osm(oi(r) — 02(r)) 2 r"^ 

— $12 cos2£:(r)} sin(6r + 5i(r)) sin(6r + 52('")) (5-27) 



Similar coupled equations are derived by for coupled S— wave equations. Since our 
potentials include the angular momentum barrier terms we use simple trigonometric func- 
tions in place of spherical Bessel and Hankel functions. This requires a modification of the 
boundary conditions just as in the uncoupled case. For this end we find it most convenient 
to rewrite the above three equations in the matrix form 

= -^[sin2(6r)$L + sin(6r) cos(6r)($LrL + Tl$l) + cos"^ {br)TL(l> lTl] (5.28) 

in which the matrix has eigenvalues of tan^i and tan52. The actual phase shifts are 

61 = 5i(r — > 00) 

62 = 62{r 00) 

e = e{r~^ 00), (5.29) 
The first boundary conditions on the above equations is 



Tl(0) = 0. 
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(5.30) 



The further numerical boundary conditions which we need for e{h), Si{h) and S2{h) are 
from (for small h) 

TL{h) = hT'^iO). (5.31) 

At small r, we can approximate our for coupled S and D states in terms of their small 
r behavior. We find 

- 4 f I ■ (5.32) 



Substitute Eg. ([Oil) and Eq.dO^) into Eq.dOg) and we find 



where 



TM = hmO) = bh I " ^ 1 (5.33) 

13 -1 + 7 



a = —7]. 



(3 = -vo 



7 = -T^- (5.34) 



45 

Then we can find e{h), tan(5_(/i) and ta.n 6^{h) by diagonalizing the matrix T^lh). The 
matrix diagonalizing T^ih) is 



cos e — sm e 
sin e cos £ 



This leads to the initial conditions 



2 

tan(2e) 



3^ 

2 , V+- 



2 2 r/-(_ 
tan5_(/i) = Til = hhi—r]- cos^e r^o cose sine — ( — I ) sin^e) 

3 3 45 



2. , 2 ,2 , r/+. 



tan (5+(/i) = T22 = 6/i(— sin eH — r^o cos e sin e — ( — I ) cos e) (5.35) 

3 3 45 



and from these initial conditions we can then integrate our equations ( p.25[ - p?2^ ) for the 



coupled system. (Note how these reduce to the uncoupled initial condition Eq.( p.l4D with 
no coupling.) 
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VI. PHASE SHIFT CALCULATIONS 



It is our aim to determine if an adequate description of the nucleon-nucleon phase shifts 
can be obtained by the use of the CTBDE to incorporate the meson exchange model. In 
contrast to the relativistic equations used in other approaches |T^, the 



CTBDE can be exactly reduced to a local Schrodinger-like form. This allows us to gain 
additional physical insight into the nucleon-nucleon interactions. We test our two models to 
find which one give us the best fit to the experimental phase shift data in nucleon-nucleon 
scattering. These two models are among many which have been tested. 

The data set |^ which we used in our test consist of pp and np nucleon-nucleon scattering 
phase shift data up to T^ab = 350 MeV published in physics journals between 1955 and 1992. 
In our fits, we use experimental phase shift data for NN scattering in the singlet states ^5*0, 
^Pi, and triplet states ^Pq, ^Pi, ^Si, ^Di. We use our parameter fit results from np 
scattering to predicate the result in pp scattering. (The variable phase method for potentials 
including the Coulomb potential is reviewed in Appendix D.) Thus we did not put the pp 
scattering data of singlet states ^5*0, and triplet states ^Pq? ^Pi iiito our fits. (There is no 
pp scattering in ^Pi, '^Si and ^Di states because of the consideration of the Pauli principle). 

We use 7 angular momentum states in our fit. There are 11 data points for every angular 
momentum state, in the energy range from 1 to 350 MeV, so the total number of data points 
in our fits is 77. To determine the free coupling constant (and the sigma mass irifj) in our 
potentials, we have to perform a best fit to the experimentally measured phase shift data. 
The coupling constant are generally searched by minimizing the quantity ■ The definition 
of our is 



s.th _ srexp 



where the 5f^ is theoretical phase shifts, the 5f^^ is experimental phase shifts and we let 
A(5j = 1 degree. (Our model at this stage is too simplified to perform a fit that involves the 
actual experimental errors). 

We have tried several methods to minimize our x^: the gradient method, grid method 
and Monte Carlo simulations. Our x^ drops very quickly at the beginning if we search by the 
gradient method, then it always hits some local minima and cannot jump out. Obviously, 
the grid method should lead us to the global minimum. The problem is that if we want to 
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find the best fit parameters we must let the mesh size be small. But then the calculation 
time becomes unbearably long. On the other hand if we choose a larger mesh, we will miss 
the parameters which we are looking for. 

We found that the Monte Carlo method can solve above dilemma. We set a reasonable 
range for all the parameters which we want to fit and generate all our fitting parameters 
randomly. Initially, the calculation time is also very long for this method, but it can leads 
us to a rough area where our fitting parameters are located in. Then we shrink the range 
for all our fitting parameters and do our calculation again (or use the gradient method in 
tandem), our calculation time then being greatly reduced. By repeating several time in the 
same way, we can finally find the parameters. 

To expedite our calculations further, we put restrictions on ^5*0 and '^5'i states. After 
every set of parameters is generated randomly, we first test it on the ^5*0 state at 1 MeV. 
For ^Sq state, if 

I <5f |>0.2 I C I (6.2) 

we let the computer jump out of this loop and generate another set of parameters and test it 
again until a set of parameters passes this restriction. Then we test it on the ^Si states at 1 
MeV with the same restriction. We only calculate at higher energy if a set of parameters 
pass these two restrictions. Our code can run at least 50 times faster by this two restrictions. 
After we shrink our parameter ranges 2 or 3 times, all of our parameters are confined in a 
small region. At this time, we may change our restriction to 

I 5f l> 0.15 I I (6.3) 

and put restriction on ^Pi states or any other states to let our code run more efficiently. 

Using this method we tried several different models to fit the phase shift experimental 
data of seven different angular momentum states including the sing let states ^S'o, ^Pi, ^£'2 
and triplet states ^Pq, ^Pi, ^Si, ^Di. Two models which we discussed above can give us a 
fairly good fit to the experimental data. The parameters which we obtained for model 1 are 
listed in table |I|, and for model 2 are listed in table |ITT| . For the features of mesons in table 
m and m, please refer to table | and Eq.( [4.4|) to Eq. (|4.6|) . The sigma mass is in MeV while 
the structure parameter tq is in inverse MeV. 
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TABLE II: Parameters From Fitting Experimental Data(Model 1). 





V 


r]' 


cr 


P 


UJ 


TT 


ao 


y 


2.25 


4.80 


47.9 


11.6 


16.5 


13.3 


0.13 


ro(xlO-3) 


2.843 


2.843 


2.843 


2.843 


2.843 


0.645 


2.843 






/o 


P' 


cu' 


<P' 


m„ 




g' 


5.64 


19.9 


0.34 


20.6 


3.10 


724.1 




ro(xlO-3) 


2.843 


2.843 


2.843 


2.843 


2.843 







TABLE III: Parameters Prom Fitting Experimental Data(Model 2). 





V 




o 


P 


LJ 


TT 


ao 




0.88 


1.70 


54.7 


2.58 


18.3 


13.6 


10.5 




1.336 


1.264 


3.180 


6.640 


2.627 


1.717 


9.282 







/o 


P' 




<\>' 


rria 




g' 


9.12 


33.5 


5.11 


28.6 


12.1 


694.3 




ro(xlO-3) 


11.45 


4.447 


6.640 


2.627 


11.45 







A. Model 1 

The theoretical phase shifts which we calculated by using the parameters for model 1 
and the experimental phase shifts for all the seven states are listed in table We use 
parameters given above to predict the phase shift of scattering. Our prediction for the 
four -p-p scattering states which include singlet states ^S'o, ^-D2 and triplet states '^Po? are 
listed in table 0. 
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TABLE IV: np Scattering Phase Shift Of ^^o, ^Pi, ^£>2, ^-Po, ^Pi, ^Si And ^Di States(Model 1). 



Energy 


^^0 


'Pi 




3 p 

Pq 


(MeV) 


Exp. The. 


Exp. The. 


Exp. The. 


Exp. The. 


1 


62.07 59.96 


-0.187 -0.359 


0.00 0.00 


0.18 0.00 


5 


63.63 63.48 


-1.487 -1.169 


0.04 0.00 


1.63 1.55 


10 


59.96 60.40 


-3.039 -2.870 


0.16 0.05 


3.65 3.57 


25 


50.90 51.95 


-6.311 -6.641 


0.68 0.52 


8.13 8.72 


50 


40.54 41.65 


-9.670 -10.23 


1.73 1.13 


10.70 11.62 


100 


26.78 26.64 


-14.52 -13.49 


3.90 2.00 


8.460 10.17 


150 


16.94 15.18 


-18.65 -15.26 


5.79 2.51 


3.690 5.688 


200 


8.940 5.615 


-22.18 -16.49 


7.29 2.91 


-1.44 0.66 


250 


1.960 -2.719 


-25.13 -17.60 


8.53 3.11 


-6.51 -4.38 


300 


-4.460 -10.16 


-27.58 -18.63 


9.69 3.55 


-11.47 -9.206 


350 


-10.59 -16.94 


-29.66 -19.68 


10.96 3.311 


-16.39 -13.81 


Energy 


'Pi 






e 


(MeV) 


Exp. The. 


Exp. The. 


Exp. The. 


Exp. The. 


1 


-0.11 -0.33 


147.747 142.692 


-0.005 0.719 


0.105 0.287 


5 


-0.94 -0.88 


118.178 112.670 


-0.183 -0.176 


0.672 1.224 


10 


-2.06 -2.26 


102.611 98.215 


-0.677 -0.256 


1.159 1.951 


25 


-4.88 -5.70 


80.63 78.38 


-2.799 -2.910 


1.793 2.587 


50 


-8.25 -10.18 


62.77 62.00 


-6.433 -6.947 


2.109 2.495 


100 


-13.24 -16.66 


43.23 43.18 


-12.23 -13.94 


2.420 3.013 


150 


-17.46 -22.12 


30.72 30.64 


-16.48 -19.35 


2.750 3.562 


200 


-21.30 -26.98 


21.22 20.95 


-19.71 -23.78 


3.130 4.489 


250 


-24.84 -31.46 


13.39 12.95 


-22.21 -27.62 


3.560 5.682 


300 


-28.07 -35.67 


6.600 6.127 


-24.14 -31.01 


4.030 6.982 


350 


-30.97 -39.58 


0.502 0.171 


-25.57 -34.15 


4.570 8.536 
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TABLE V: pp Scattering Phase Shift Of ^Sq, ^D2, ^Pq And ^Pi States(Model 1). 



Energy 


^5*0 






'Pi 


MeV 


Exp. The. 


Exp. The. 


Exp. The. 


Exp. The. 


1 


32.68 51.95 


0.001 -0.091 


0.134 0.381 


-0.081 -1.215 


5 


54.83 55.47 


0.043 -0.183 


1.582 0.954 


-0.902 -2.536 


10 


55.22 54.45 


0.165 -0.270 


3.729 1.773 


-2.060 -3.864 


25 


48.67 47.64 


0.696 -0.441 


8.575 5.422 


-4.932 -7.932 


50 


38.90 37.77 


1.711 -0.504 


11.47 9.766 


-8.317-13.15 


100 


24.97 23.63 


3.790 0.511 


9.450 7.862 


-13.26 -18.45 


150 


14.75 12.37 


5.606 1.141 


4.740 3.812 


-17.43 -24.42 


200 


6.550 3.024 


7.058 2.407 


-0.370 -1.178 


-21.25 -28.50 


250 


-0.31 -5.15 


8.270 2.994 


-5.430 -6.193 


-24.77 -33.26 


300 


-6.15 -12.55 


9.420 3.136 


-10.39 -10.98 


-27.99 -37.63 


350 


-11.13 -19.27 


10.69 2.902 


-15.30 -15.42 


-30.89 -41.13 



The resuhs for np scattering are also presented from figure |l| to figure ^ and for pp 



scattering from figure a to figure 11 
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FIG. 4: np Scattering Phase Shift for ^Pq State (Model 1) 
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FIG. 6: np Scattering Phase Shift for State (Model 1) 
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FIG. 9: np Scattering Phase Shift for State (Model 1) 
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FIG. 10: pp Scattering Phase Shift for ^Pq State (Model 1) 
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B. Model 2 



The theoretical phase shifts which we calculated by using the parameters for model 2 and 
the experimental phase shifts for all the seven states are listed in table [V^. We also use the 
parameters for model 2 to predict the phase shift of pp scattering. 
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TABLE VI: np Scattering Phase Shift Of ^^o, ^Pi, ^£>2, ^-Po, ^Pi, ^Si And ^Di States(Model 2). 



Energy 


^^0 


'Pi 




3 p 

Pq 


(MeV) 


Exp. The. 


Exp. The. 


Exp. The. 


Exp. The. 


1 


62.07 60.60 


-0.187 -0.358 


0.00 0.02 


0.18 0.00 


5 


63.63 63.50 


-1.487 -1.163 


0.04 0.15 


1.63 1.61 


10 


59.96 60.20 


-3.039 -2.857 


0.16 0.39 


3.65 3.74 


25 


50.90 51.44 


-6.311 -6.629 


0.68 0.40 


8.13 9.28 


50 


40.54 40.91 


-9.670 -10.36 


1.73 1.37 


10.70 12.69 


100 


26.78 25.86 


-14.52 -14.44 


3.90 2.42 


8.460 11.74 


150 


16.94 14.62 


-18.65 -17.55 


5.79 3.62 


3.690 7.399 


200 


8.940 5.435 


-22.18 -20.37 


7.29 4.55 


-1.44 2.36 


250 


1.960 -2.428 


-25.13 -23.15 


8.53 5.24 


-6.51 -2.78 


300 


-4.460 -9.330 


-27.58 -25.87 


9.69 5.34 


-11.47 -7.746 


350 


-10.59 -15.52 


-29.66 -28.54 


10.96 5.30 


-16.39 -12.52 


Energy 


'Pi 






e 


(MeV) 


Exp. The. 


Exp. The. 


Exp. The. 


Exp. The. 


1 


-0.11 -0.32 


147.747 144.797 


-0.005 0.719 


0.105 0.264 


5 


-0.94 -0.81 


118.178 115.232 


-0.183 -0.172 


0.672 1.106 


10 


-2.06 -2.08 


102.611 100.668 


-0.677 -0.239 


1.159 1.723 


25 


-4.88 -5.07 


80.63 80.66 


-2.799 -2.834 


1.793 2.099 


50 


-8.25 -8.68 


62.77 64.30 


-6.433 -6.798 


2.109 1.708 


100 


-13.24 -13.55 


43.23 45.68 


-12.23 -13.77 


2.420 1.663 


150 


-17.46 -17.74 


30.72 33.35 


-16.48 -19.34 


2.750 1.541 


200 


-21.30 -21.67 


21.22 23.80 


-19.71 -24.11 


3.130 1.648 


250 


-24.84 -25.47 


13.39 15.90 


-22.21 -28.38 


3.560 1.834 


300 


-28.07 -29.14 


6.600 9.099 


-24.14 -32.29 


4.030 1.965 


350 


-30.97 -32.67 


0.502 3.095 


-25.57 -36.01 


4.570 2.147 
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TABLE VII: pp Scattering Phase Shift Of ^Sq, ^Ds, ^^o And ^Pi States(Model 2). 



Energy 








'Pi 


MeV 


Exp. The. 


Exp. The. 


Exp. The. 


Exp. The. 


1 


32.68 52.40 


0.001 -0.116 


0.134 0.417 


-0.081 -1.172 


5 


54.83 55.48 


0.043 -0.232 


1.582 1.042 


-0.902 -2.434 


10 


55.22 54.24 


0.165 -0.327 


3.729 1.934 


-2.060 -3.682 


25 


48.67 47.13 


0.696 -0.524 


8.575 5.943 


-4.932 -7.355 


50 


38.90 37.04 


1.711 -0.505 


11.47 10.88 


-8.317 -11.57 


100 


24.97 22.85 


3.790 0.994 


9.450 9.417 


-13.26 -15.41 


150 


14.75 11.82 


5.606 2.036 


4.740 5.543 


-17.43 -19.97 


200 


6.550 2.845 


7.058 3.211 


-0.370 0.495 


-21.25 -23.23 


250 


-0.31 -4.86 


8.270 3.648 


-5.430 -4.589 


-24.77 -27.28 


300 


-6.15 -11.72 


9.420 3.956 


-10.39 -9.516 


-27.99 -31.05 


350 


-11.13 -17.85 


10.69 4.014 


-15.30 -14.13 


-30.89 -34.22 



The prediction for the four pp scattering states are Hsted in table VD. The results 



of model 2 for np scattering are given in figure 12 to figure ^ and for pp scattering are 



from figure |T9| to figure |22|. Our results show for this model an improvement over those of 
model 1 especially for the the singlet P and D states. However there is still much to be 
desired in the fit. One possible cause of this problem is that we did not include tensor and 
pseudovector interactions in our covariant potentials, limiting ourselves to scalar, vector and 
pseudoscalar. Another may be the ignoring of the pseudovector coupling of the pseudoscalar 
mesons to the nucleon. Our results in pp scattering show that if we obtain a good fit in 
np scattering our predicted results in pp scattering will also be good. This means that it 
is unnecessary to include pp scattering in the our fit, we may use the parameters obtained 
in np scattering to predict the results in pp scattering. Overall our results are promising 
and indicate that the two-body Dirac equations of constraint dynamics together with the 
meson exchange model are suitable to construct semi-phenomenological potential models for 
nucleon-nucleon scattering. 
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FIG. 12: np Scattering Phase Shift of ^Sq State (Model 2) 



VII. CONCLUSION 

The two-body Dirac equations of constraint dynamics constitutes the first fully covariant 
treatment of the relativistic two-body problem that 

a) includes constituent spin, b) regulates the relative time in a covariant manner, c) 
provides an exact reduction to 4 decoupled 4-component wave equations, d) includes non- 
perturbative recoil effects in a natural way that eliminates the need for singularity-softening 
parameters or finite particle size in semiphenomenological applications to QCD, e) is canon- 
ically equivalent in the semi-relativistic approximation to the Fermi-Breit approximation 
to the Bethe-Salpeter equation, f) unlike the Bethe-Salpeter equation and most other rela- 
tivisitc approaches has a local momentum structure as simple as that of the nonrelativistic 
Schrodinger equation, g) is well-defined for zero-mass constituents (hence, permits investi- 
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FIG. 13: np Scattering Phase Shift of ^Pi State (Model 2) 



gation of the chiral symmetry limit) h) possesses spin structure that yields an exact solution 
for singlet positronium, i) has static limits that are relativistic, reducing to the ordinary 
single-particle Dirac equation in the limit that either particle becomes infinitely heavy, j) 
possesses a great variety of equivalent forms that are rearrangements of its two coupled Dirac 
equations (hence is directly related to many previously-known quantum descriptions of the 
relativistic two-body system). These structures play an essential role in the success of this 
approach to both QCD and QED bound states. What is noteworthy in the latter applhca- 
tion is that one need only identify the non-relativistic parts, i.e. the lowest order forms of A 
and S. The spin-dependent and covariant structure of the two-body Dirac formalism then 
automatically stamps out the correct semirelativistic spin dependent and spin independent 
corrections and provides well defined higher order relativistic corrections as well. In addi- 
tion the constraint formalism, although rooted in classical mechanics, has close connections 
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FIG. 14: np Scattering Phase Shift of State (Model 2) 



to the Bethe-Salpeter equation of quantum field theory [^] and with Wigner's formulation 
of relativistic quantum mechanics as a symmetry of quantum theory [ ^ . 

In this paper we have shown that these two-body Dirac equations may provide a rea- 
sonable account of the nucleon-nucleon scattering data when combined with the meson 
exchange model. What makes this result important is that it is accomplished with a lo- 
cal and covariant formulation of the two-body problem. What makes this unique is that 
this approach has been thoroughly tested in a nonperturbative context for both QED and 
QCD bound states. It is not a given that success in one or even both areas would imply 
that the formalism would do well in another. In particular, the fits could have easily been 
disastrous given that the minimal coupling idea we have used (based in part on the earlier 
work on the quasipotential approach of Todorov). The reason for some doubt is that these 
minimal coupling forms (generalized to the scalar interactions as well as the vector) lead 
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FIG. 15: np Scattering Phase Shift of ^Pq State (Model 2) 

to the scalar and vector potentials appearing squared. Because of the size of the coupling 
constants, the deviation from the standard effective potentials could have been considerable 
in all cases. There are other nonperturbative structures that appear in the Pauli reduction 
of our equations to Schrodinger-like form (typical of what appears in the Pauli reduction of 
the one-body Dirac equation) that could also have prevented any reasonable results. So the 
general agreement wc obtained with the data is very encouraging that this approach could 
be extended to include more general interactions. 

An important step in our reduction was that we put the equation in a form for which we 
can apply the techniques which have been already developed for the Schrodinger-like system 
in nonrelativistic quantum mechanics. This required that we get rid of first derivative 
terms. For the uncoupled states, it is pretty straightforward. For the coupled states we 
used a different spin-matirx approach that works for both the uncoupled and coupled states 
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FIG. 16: np Scattering Phase Shift of ^Pi State (Model 2) 



simultaneously. 

We then tested several models by using the variable phase methods. We found it 
most convenient to put all the angular momentum barrier terms in the potentials, and 
change all the phase shift equations to the form of S state-like phase shift equations (see 
Eqs. ([5lI| , p:25| , p:26ip:27D ) . 

After several models and several methods to minimize our tested, we found two models 
which can lead us to a fairly good fit to the experimental phase shift data. 

The most important equation used in our phase shift analysis for nucleon-nucleon scat- 
tering is Eq.(g:32D . It is a coupled Schrddinger-like equation derived from two-body Dirac 
equations with no approximations. All of our radial wave equations for any specific angular 
momentum state are obtained from this equation. 

We use nine mesons in our fit. We summarize the meson-nucleon interactions we used by 
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FIG. 17: np Scattering Phase Shift of State (Model 2) 

writing the quantum field theory Lagrange function for their effective interactions 

Li = Qa^ipo- + gf^i/Ji/jfo + gaoi^r^p ■ ao 

+ gAl^Til) ■ TT + gr^ip^^i/^rj + grf'tp-f^i/jri' (7.1) 

where ip represent the nucleon field, a, /o, ... represent the meson fields. 

Several models have been tested by using the variable phase methods, two models can 
lead us to a fairly good fit to the experimental phase shift data. We use the parameters which 
gives good fits to the np scattering data to predict the phase shifts for the pp scattering. 
These lead to a good prediction for the pp scattering based on the parameters we obtained 
(with noted exceptions). This means that our work has shown a promising result. The 
following are some suggestions to improve our work in the future. 
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FIG. 18: np Scattering Phase Shift of ^Di State (Model 2) 

VIII. SUGGESTIONS FOR FUTURE WORK 
A. Other Model Tests 

More model testing is absolutely necessary in the future. By model we mean the way we 
place the perturbative interactions that arise from Eg. (|7. 1| ) into the nonperturbative forms 
we need for L, C and Q. During our fits, we found that our final result are sensitive to the 
model we chose ranging from very bad fits to the fits presented here. Changing the way to 
modify the interactions and the way of mesons enter into the two-body Dirac equations may 
provide a new opportunity to improve our fit. 
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FIG. 19: pp Scattering Phase Shift of State (Model 2) 



B. Including World Tensor Interactions 

We have included just scalar, pseudoscalar and vector interactions in our potentials 
through the invariant forms like L, C and Q. Treating two-body Dirac equations with 
tensor interactions of the vector meson may improve our fit. These tensor interacting were 
discussed earlier (see Eq.( [1.16D ) and correspond to non- minimal coupling of spin one-half 
particle not present in QED but which can not be ruled out in massive vector meson-nucleon 
interactions. The corresponding field theory interaction is 

AL, = g';;^a^''T,ljp^, + g'^Pa^^'^iJuj^, + g'^i^a^^'^Jct^^,, (8.1) 
and would correspond to relaxing the free field equation assumption made in Eq. ( [4.17|) ). 
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FIG. 20: pp Scattering Phase Shift of State (Model 2) 



C. Include Pseudovector Interactions 

Another option is to allow the pseudoscalar mesons(7r, rj, and rj') to interact with the 
nucleon not only by the pseudoscalar interaction(as in Eg. ( |7.1| )) but also by the way of the 
pseudovector interactions as below 

ALi = g'JPj^j^'Tijd.n + g'^i^Y^^'^jd^r, + ^;,^7^75^9^r/' (8.2) 

D. Include Full Massive Spin-One Propagator 

We have ignored a portion of the massive spin-one propagator in our fit which is zero 
for particles on the mass shell. To include this portion of massive spin-one propagator we 
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FIG. 21: pp Scattering Phase Shift of ^Pq State (Model 2) 



would have to change the vector propagator as below 



(8.3) 



+ — ie + — ie 

Among all the four suggestions, the first one would be technically easiest once we finds 
models more general than the two we have presented here . The last three suggestions 
would involve corresponding additions to the interaction that appear in the two-body Dirac 
equations. Because the above interactions all involve derivative couplings we will have to 
examine the CTBDE for the corresponding invariant A's. These would include not only 
the eight invariants listed earlier (see Eg. ( |2.46| ) to Eq.( p^.56D ) but also four additional ones 



corresponding to A = R9i ■ £±02 ■ xj_, 2S0^i0^26i ■ £±02 ■ x± , 2T6i ■ P9i ■ POi ■ £±62 ■ x± 
and 4f/6'5i6'520i ■ P9i ■ POi ■ xj_92 ■ x_i_. The four functions R, S,T, U are each functions of 
x± and they represent space-like interactions paralleling those corresponding to Q, I, Y, and 
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FIG. 22: pp Scattering Phase Shift of ^Pi State (Model 2) 



respectively given earlier. To include all 12 covariant matrix interactions will involve 
a significant modification of our basic equation Eq. (|3.32|) as well as the two-body Dirac 
equations given in Eqs( |2.58| , |2.59|) . 



E. Extentions to the N-Body Problem 



Can the constraint formalism be extended to iV-bodies? There is no solution to the 
compatibility condition 

[n,,n,]\^) = 0; t,j = l,..,N (8.4) 

of generalized mass-shell constraints (or their Dirac counterparts) that has the simplicity 
of the "third law" and tranversality conditions given in ( gl4D and (^JEj) . The difficulty 
involves satisfying Eq. (|8.4|) and cluster separability (needed to describe scattering states) 
at the same time. Rohrlich has shown that this necessarily involves the introduction of 
A^— body forces [^]. If one is willing to limit A^— body considerations to bound states (so 
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that cluster considerations are not important) then Ref. [|53| provides a constraint formahsm 



in which a single dynamical wave equation (as in the two-body case) determines the bound 
state energies. Ref. (and references contained therein) provides an A^— body constraint 
formalism that involves particles and fields leading in the end to directly inteacting particles 
by elimination of the field degrees of freedom by second class contraints. 

APPENDIX A: PAULI-FORM OF TWO-BODY DIRAC EQUATIONS 



We rewrite Eqs.( p.6l] , p.62| ) by multiplying the first by v^i/?i and the second by y/2i(32 
yielding B 



[TiiW + t/i(/3i/52)75i752]|^) = (^1 + MiA)7, 
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-[T2(A/?2) + ?72(A/32)751752]|^) = (^2 + (A.l) 

in which the kinetic and recoil terms are 

Ti(/3i/?2) = exp(6;)[Si ■ p - ^/3i/32(S2 ■ d{-C + GP.P^^, ■ S2)] 



T2{W = exp{g)[^2 ■ P - '-PM^i ■ d{-C + g(3,(32^, ■ S2)] (A.2) 
f/i(A/52) = exp(g)[-^A/?25]2 ■ a(J/?i/?2 - L)] 



U2{PiP2) = exp(6;)[-^/5i/32Si ■ d{jPiP2 - L)] (A.3) 



while the timelike and scalar potentials Ei, Mi are given above in Eqs.(|2.63|) and 



The final result of the matrix multiplication in Eqs. (|A.l| ) is a set of eight simultaneous 
equations for the Dirac spinors \ip)2, |V')3; |V')4- arbitrary frame, the result of the 

matrix calculation produces the eight simultaneous equations {<7i\ip) I V") 1,2,3,4) HI- One 

then reduces the eight equations to a second order Schrodinger-like equation by a process of 
substitution and elimination using the combinations of the four Dirac-spinors given below 
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|0±)^|^)l±|^)4, 



\X±) = 1^)2 ± (A.4) 
We display all the general spin dependent structures in $(r, p, cri, 0-2, ty) explicitly, very 
similar to what appears in nonrelativistic formalisms such as seen in the Hamada- Johnson 
and Yale group models (as well as the nonrelativistic limit of Gross's equation). We do 
this by expressing it explicitly in terms of its matrix (cri,cr2), and operator p structure 
in the CM system {P = (1,0) ). We are working in the CM frame (i.e. x± = (r,0)), so 
all the interaction functions {L{x±), J{x±),C{x±),Q{x±)) are functions of r = \/x^ = |r|, 
F = F{r) 

Ref. PI finds the reduction 

h[Ei[(Ti ■ p - icT2 ■ (d + kcTi ■ cr2)]hFi[cri ■ p - i(T2-{z + kcri ■ cr2)]|0+) 
+h[Mi[(Ti ■ p - i(T2-{o + kcTi ■ (T2)]hF3[cri ■ p - icr2-(z + kcri • cr2)]|0+) 
-h[Ei[cri ■ p - icr2-(d + kcri ■ (T2)]hF2[(T2 ■ p - i(Ti-{z + kcri ■ cr2)]|0+) 
+h[Mi[(7i ■ p - i(T2-{o + kcTi ■ (T2)]hF4^[(T2 ■ p - icri-{z + kcTi ■ 0-2)] 10+) 



= fi>+). (A.5) 

in which 

B"^ = Ei^ - Mi^ = E2^ - M2^ 

= b'^iw) + (ej + e^)sinh2( J) + 2eie2 sinh( J) cosh( J) (A.6) 
— {ml + m2)sinh^(L) — 2mim2 sinh(L) cosh(L). 

and 

h = exp(^), 
k^lviog(/i), 
z= ^V(-C + J-L) 
d=^\/{C + J + L) 

o = \v{C - J - L). (A.7) 



71 



with 



M2 
V 

^ V 
F3 = — 



(A.8) 



V = EiM2 + E2Mi. (A.9) 



Eq.( [A.5| ) is a second-order Schrodinger-like eigenvalue equation for the newly defined wave- 



function 10+) in the form. 

{pl + <l>{r,(T,,(T2,w))\<P+) =b\w)\<l)+). (A.IO) 
Eq.( |3.5|) for B"^ provide us with the primary spin independent part of $, the quasipotential. 



Note that in the CM system p\ = p^, a = (0,cr).For future reference we will refer to the 
four sets of terms on the left hand side as the Eq. (|A.5|) (a),(b),(c),(d) term. 

Now we proceed with a different derivation than Long and Crater's derivation The 
aim is to produce a Schrodinger like form like in Eq.( |A.10|) involving the Pauli matrices for 
both particles. 



Substitute d, h, Fi, z, k's expressions to (a) term of Eq. (|A.5|) , we obtain 



(a) term = exp(^)Ei{[<Ti ■ p - ^0-2 ■ V(C + J + L) - ^VG ■ (cri + ien x cts) 



X exp(e?)^[o-i ■ p - ^0-2 ■ V(-C + J - L) - ■ (cti + icr^ x 0-2)]} (A.ll) 

working out the commutation relation of cti ■ p in above expression, we can find the (a) term 
is 

(a) term= exp(^)_E'iX 

M2 ^ i i 
{exp(^) — [p --cr2-V(-C+J-L)(cri-p)--V^-[(p-fi(criXp)-(cri-cr2)p+cri(cr2-p)-i(cr2Xp 

+ ■ 9[exp(^)^[o-i . p - lcr2 ■ V(-C + J - L) - ^V^ ■ (cti + icr^ x 0-2)]]- 
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^[<T2-V(C+J+L)+Ve?-(<Ti+i<TiXo-2)]exp(^)^[CTrp-^o-2-V(-C+J-L)-^V^-((Ti+itriXfr2)]^ 

Likewise we can the find (b),(c),(d) terms, 
(b) term= exp(^)MiX 



E ^ i i 
{exp(6?)^[p --(T2-V(-C+J-L)((Trp)--V^-[(p+i(o-iXp)-((Tr(T2)p+<Ti((T2-p)-i(<T2Xp)]] 



■ d[exp{g)^[ai ■ p - ^0-2 ■ V(-C + J - L) - ^VQ ■ {en + iai x 0-2)]]- 



^[o-2-V(C-J-L)+V^-(o-i+io-iXo-2)]exp(^):^[o-rp-^o-2-V(-C+J-L)-^V^-(ai+io-i^ 



(c) term= — exp(^)ii^iX 



{exp(^)^[(<T2-p)(<Ti-p)-^<TrV(-C+J-L)(o-i-p)-^V^-[((T2(<Trp)-(<Tr<T2)p+o-i(<T2-p) 



+^(0-2 X p)]] + \(Ti ■ 9[exp(^)^[<r2 • p - Jo-i • V(-C + J - L) - ^VG ■ ((T2 + i<T2 x (Ti)]]- 



i Ml i i 

-[<r2-V(C+J+L)+V^-(<Ti+i<TiX<72)]exp(^) — [(T2-p--<TrV(-C+J-L)--Ve?-(o-2+io-2X<Ti)]} 

(d) term= exp(^)MiX 



{exp(^)^[(o-2-p)(o-rp)-^<TrV(-C+J-L)(o-rp)-^V^-[(o-2(<Trp)-(o-rO-2)p+o-i(o-2-p) 



+i(<T2 X p)]] + -o-i • a[exp(^?)|i[o-2 • p - ^o-i • V(-C + J - L) - ^V^ • (0-2 + i(72 X o-i)]]- 



^[a2-V{C-J-L)+Vgiai+i(riXa2)]exp{g)^[a2-p-^(TvV{-C+J-L)-^Vg-{(T2+i(T2Xai)]} 
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We simplify the above expressions by using identities involving cri and (J2 and group 
above equations by the term , Darwin term (r ■ p), spin-orbit angular momentum term 
L ■ {(Ti + 0-2), spin-orbit angular momentum difference term L • (cti — (T2), spin-spin term 
(cTi ■ CT2), tensor term {cri ■ r){cr2 ■ f), additional spin dependent terms L ■ {(Ti x 0-2) and 
(cTi • r)(cr2 • p) + {(T2 ■ r){(Ti ■ p), and spin independent terms. Collecting all terms for the 
(a) + (b) + (c) + (d) terms our Eq.([0|) becomes Eq.Qg^). 

APPENDIX B: RADIAL EQUATIONS 

The following are radial eigenvalue equations corresponding to Eq.(|3.4|) after getting rid 
of the first derivative terms for singlet states ^Sq, ^Pi, ^D2{ a general singlet ^Jj), triplet 
states ^Pi( a general let ^Jj), a general s = 1, j = / + 1 ( ^Po,^Si states ), and a general 
s = 1, j = / + 1 (^Di state). 

^Sq, ^Pi, ( a general singlet ^Jj) L ■ (cri + (T2) = 0, (Ti ■ (T2 = —3, cti ■ r(T2 ■ r = —1. 



{-^ + ^^^^^ + ^4 + h" + ^-!^ + --3k-j-g'h'-h"~ — + m}v = B'exp{-2g)v (B.l) 
ar r'^ A 1 r r 

^Pi{ a general triplet ^Jj) L ■ {cri + (T2) = —2, cri ■ (T2 = 1, cri ■ i<^2 ■ ? =1- 



{-^ + ^^-^t + ^ + + C + ^ + ^ + 9'h' + h" + m}v = B"" eM-'^Q)v, (B.2) 

s = 1, j = / + 1 ( states ) 

L ■ (cTi + 0-2) = 2(j — 1), cri ■ a2 = 1, cri ■ rcr2 ■ r =2fpi (diagonal term), and cri ■ fcr2 • 
f = ^^j!^i^ - (off diagonal term). 



d'^ 7(7 — 1) Sg'sinh^/i , , cosh /i sinh /i 6sinh^/i gf' cosh /i sinh 
{-T^ + o +— + 5 



, , sinh^ h cosh h sinh qf'^ , q" q' 
-2h' + 2 + ^ + h'^ + ^ + - + k 



,,q' fif'sinh^/i sinh^ /i ,, cosh sinh /i , 
-F2 J - 1 f + ^ 2 — + 2h' ] 

^^y fy ry^ ry 
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2(7 — 1) r^, /Sinh h ^ cosh /i sinh /i h' g'cosh/isinh/i, 

+ y^/ 2^ 2 — - — + - + ^ 

1 . 3g' cosh /i sinh h q' sinh^ /i , , cosh h sinh /i , , sinh^ h 

+- [— 2h' + 6h' 

2j + 1 r r r r 

cosh /i sinh /i ^sinh^/i ,,, , ,, 2/i', , 
-6 h 2 — + n + + h" H + m}u+ 

2\/ jij + 1) fSg' cosh /i sinh /i g'sinh^/i ^, .cosh /i sinh /i ^,,smh^h 

H — ^^^^-^^^ -{— 2h' h 6/i' 

2 j + 1 r r r r 

cosh /i sinh /i sinh^ /i , . 2/i' 

-6 . + 2 — + n + ^'/i' + /i" + — 



+2(i - _ 2cosh(^^)sinh(/.) ^ ^ ^ ,-cosh(/.^)sinh(/.)j^^_ ^ ^2 ^^p(_2^)^^^ 

(B.3) 

s = 1, j = / - 1 ( 3Po,32p^ g^^^gg ) 

L • (<Ti + (T2) = -2{j + 2), cTi • (T2 = 1, fTi • f<T2 • f = -^^(diagonal term), and <Ti • f<T2 • 

f = ^^^^^j'" - (off diagonal term). 



d'^ (j + l + 2) 3g' siidr /?. .cosh /; sink /i 6 sinh^ /i g' cosh /; smh /i 
{-^ + + 6^ ^ 7- r 



,,sinh^/i „ cosh /i sinh /i g'^ g" g' 

-2/i' + 2 + ^ + h'^ + ^ + - + k 

r r"' 4 2 r 



„, „vr9' g'sinh^/i sinh^ /i . cosh /i sinh /i, 

+20 + 2 f + ^ 2 ^ + 2h' 

2r r r r 



2(7 — 1) r^, ,sinh^ /i ^ cosh /i sinh /i /i' g' cosh /i sinh /i, 

+^j--r 2/i' 2 ^ — + - + ^ 

1 ,3g' cosh /i sinh /i g'sinh^/i „, .cosh /i sinh /i „,,sinh^/i 
^ ^ ^ 2/i' h 6/i' 



2j + 1 r r 

cosh /t sinh /t sinh^/i . , , ,, 2h' . 
-6 h 2 h n + c/'/i' + h" H + m}w_ 
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2 a/jXj + 1) 1 3 (?' cosh h sinh /i 



2j + l 



fif' sinh^ /i , , cosh h sinh /i , , sinh^ h 
2h' h 6/i 



cosh /i sinh /i sm\i' h ,, , , ,, 2h' 

-6 5 + 2 — + n + g'h' + h" + — 



-2(j + 2)[^^^1^^^ - 2cosh(/.)sinh(/.) ^ ^ ^-cosh(/.)sinh(/.) ^ exp(-2^)«_, 

(B.4) 

Substituting for (?', m, n, /c we obtain the radial equations and potentials $ given in 
the text. 



APPENDIX C: DERIVATION OF COUPLED PHASE SHIFT EQUATIONS 



We have found that we can use the Messiah ansatz [H6| 



for the solution of 



to yield 



u = A sin(6r + 5) 

u' = bAcos{br + 5) (C.l) 

■^ + ^L{r))u = b\{r) (C.2) 

5' = + (C.3) 



Next we see how this can be worked out in the case of coupled radial equations of the 
form 

bA""- |. (C.4) 

where $l is a two by two matrix. This equation will have solutions that are S*— wave 
dominant and wave dominant. Form them together into a 2 x 2 matrix U so that the 
above equation becomes 

-U" + ^lU = b'^U. (C.5) 
Then take its transpose and add the two. One obtains 

-{U" + U''^) + {^lU + U^^l) = b\U + U^) (C.6) 
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In analogy to the uncoupled case we assume 

U = A{r)sm{br + V{r)) (C.7) 

where 

P=5(r) + D(r) • o- (C.8) 
A=a{r) + A(r) ■ cr. 

Let i? be a matrix that diagonalizes the phase shift matrix function V{r) to the form 

U=RUR-^ = A sm{br + S + Ddg) (C.9) 

where 

A = RAR-^ (C.IO) 

Continuing the analogy we let 

U' = bAcos{br + V). (C.ll) 

Then 

RU'R-^ = bA cos{br + 6 + Dag) = R{R-^UR)'R-^ (C.12) 
= RR-^'A sin(6r + 5 + Das) + A' sin(6r + 5 + Das) 
+ A{b + 6' + D'as) cos(6r + 6 + Das) + isin(6r + 6 + Das)R!R-^ 



But RR ^' = —R'R ^ so that we obtain the condition 

A' sm{br+S+Das)+A{S'+D'as) cos{br+S+Da3) + [Asm{br+S+Da3), R'R-^]_ = 0. (C.13) 
In general we would take 

A = a + Asas + A_l ■ cr (C.14) 
= A|| + A_L ■ cr 

and decompose Eq.( p.l3 | ) and Eq. (|C.6|) into two sets of four coupled equations. 
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Give i? the following general form: 

R = exp{ie{r)o-2) exp(?7(r)(Ti) 
R~^ = exp(— r^cTi) exp{—iea2) 
R! = ie' 02 exp(i£:o"2) exp(77(Ti) + exp(i£:cr2) exp(r7o"i)?7'ai 
R'R~^ = ie'a2 + exp{iEa2)i]'ai exp{—ie(T2) 

= ie'a2 + rl co^{2e)oi + 77' sin(2e)a3 (C.15) 

We consider the case in which is a symmtric matrix and furthermore that as a result 
V = V^ . In that case our matrix is R is orthogonal (r^ = 0). 

Next we examine the three terms of Eg. ( |C.6|) . We assume that A_l is symmetric so that 

= and 

\?R(JJ + \J'^)R-^ = + Ax ■ <t) sin(5r + 6 + Da^) + sin(6r + 6 + Da3){A\\ + • (t]) 

= 62[2i||sin(6r + 5 + D(T3) + 2(iisin(6r + 5)cosD)ai] (C.16) 

and 

R{^U + U'^^'^)R-^ 

= ($11 + $_L-cr)(A|| + A_L ■ cr) sin(6r + 5 + Da^) + (Transpose) 
= + $_L ■ A_l) sin(6r + 5 + Da^) 

+ (I'liAx • cr + $_L ■ cr^ll + x A_l ■ cr) sin(6r + 5 + Da^) 
+ sin(6r + 5 + Da3){^\\A\\ + l*^ ■ A^) 

+ sin(6r + 5 + Da3){k^ ■ cr<l|| + A\\^^-a + zI>_l x A_l ■ cr) (C.17) 

The term $x x A^ ■ cr is zero since A2 = Q = ^2- The second derivative term is 

R{U" + f/"'^)i?^ 
= R{R^U'R)'R-^ + (Transpose) 
= bR{R^Acos{br + 5 + Das)RyR^ + (Transpose) 
= b{ [A cos(6r + 5 + Das), R'R^] _ + A' cos(6r + 5 + Da^) 

- A{b + 6' + D'a-i) sm{br + 6 + /^crg)} + (Transpose) 

= b{ie'[{A\\ + Ax • <t) cos(6r + 6 + Das), (^2]- + (i|| + A^ ■ o") cos(6r + 5 + Z^ag) 

- (i|| + A±-a){b + 5' + D'as) sin(6r + 5 + Das)} + (Transpose) (C.18) 
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Using properties of the Pauli matrices and dividing Eg. ( p. 131) and Eg. (|C.6|) into || and 
components we obtain the four eguations 

- A\\{S' + D'a^) sin(6r + S + Da-,) 

+ ij| cos(6r + 6 + Das) - ^be'a^Ai cos(6r + 5) cosD 



i(<i>l|A|| + $iAi) sin(6r + S + Das), 



(C.19) 



Aji sin(6r + 6 + Da^) + A\i{6' + D'a^) cos(6r + 6 + Da^) - ^e'a^Ai sin(6r + 5) cosD 
= 0, (C.20) 



cos(6r + 6) cos DA[ — {6' sin(6r + 6) cos D + D' cos(6r + S) sin D)Ai 
+ 2e'{As cos(6r + 6) cos D — a sin(6r + 6) sin D) 
= - [0 sin(6r + 6) cos DAi — $3 cos(6r + 6) sin DAi 
+ (asin(6r + 6) cosD + y43(cos(6r + 6) sinD)$i], 



(C.21) 



A[sm{br + S) cosD + Ai{S' cos{br + S) cosD — D' sin{br + S) sinD) 
+ 2e'{a cos(6r + 6) sin D + A^ sin(6r + 6) cos D) 

= 0. (C.22) 
Combining Eg.( |Cl9D and Eg.( |C2(l| ) we obtain 



i|l (5' + D'as) - 2e'Ai sin D cos D 
i(<l||i|l + <liii) sin2(6r + 5 + Dcxg). 



Combining Eg. (10211) and Eg. ([02^ ) gives 



A16' csc(6r + (5) cos D — 2e' csc(6r + 6) sin D 
= -[((/)sin(6r + S) cosD — $3 cos(6r + 6) sinD)y4i 
+ (asin(6r + 6) cosD + A3(cos(6r + 6) sinD)$i]. 



(C.23) 



(C.24) 
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Rewrite the above two equations as 



i|l [{6' + D'a^) + i<l|| sin2(6r + 5 + Da^)] 
+ Ai{2£' sin D cos + ^l>i sin^(6r + 5 + Da^)) 
= 0, 



(C.25) 



-(asin(6r + 6) cosD + A3 cos(6r + 6) sinD)$i + 2e'acsc(br + 6) sinD 
b 

+ y4i[-(0sin(6r + 6) cosD — $3 cos(6r + 6) sinD) — S' esc {br + 6) cosZ^] 
= 0. 



(C.26) 



The first of these two equations is actually two equations 

a[S' + ^0(1 - cos 2(6r + S) cos(2L))) + ^$3 sin 2(6r + S) sm{2D)] 
AslD' + ^'^'3(1 - cos2(6r + 6) cos(2L))) + ^0sin2(6r + 6) sin(2D)] 
+ Ai{{e' sin2D + ^$1(1 -cos2(6r + 5) cos 2D)) 



and 



(C.27) 



a[D' + ^^-sll - cos2(6r + 6) cos(2D)) + -'^0sin2(6r + 6) sm{2D)] 
+ A3[6' + ^0(1 - cos2(6r + 6) cos(2i:))) + ^<l3sin2(6r + S) sin(2D)] 
+ Ai^^i sin2(6r + 5) sin 2D 
= 0. 



(C.28) 



So now together with Eq. (|C.26D 



a[- sin(6r + 6) cos(-D)$i + 26' csc(6r + 6) sin D] 
b 

+ Asl^ cos(6r + 6) sm{D)^i] 

+ Ai[-{(j) sm{br + 6) cosD — $3 cos(6r + 6) sinD) — 5'csc(6r + 6) cosD] 

= (C.29) 

we have three homogeneous equations in a, A3, Ai. We simplify these equations further by 
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assuming that Ai = 0. 



= a[6' + ^0(1 - cos 2(6r + 6) cos(2D)) + ^"Ig sin 2(6r + 6) sin(2D)] 
A3[D' + ^^'3(1 - cos2(6r + 6) cos(2D)) + ^0sin2(6r + S) sm{2D)] 
= (C.30) 



= a[D' + ^^'3(1 - cos2(6r + 6) cos(2L))) + ^0sin2(6r + 6) sin(2D)] 
+ A3[5' + ^0(1 - cos2(6r + 5) cos(2L))) + ^$3 sin2(6r + 5) sin(2D)] 
= (C.31) 



sin(6r + 5) cos(D)$i + 2e' csc(6r + 5) sin D] 
+ A^l^ cos(6r + 6) sin(£))<li] 

= (C.32) 



The solution we seek is 



6' + ^0(1 - cos 2(6r + 6) cos(2D)) + sin 2(6r + 6) sm(2D) = (C.33) 
2b 2b 



1 ~ 1 
D' + ^'^'3(1 - cos2(6r + 6) cos(2D)) + — 0sin2(6r + 6) sm{2D) = 0. (C.34) 



Let 



6=1(6, + 62) 

= ^(6,-62), (C.35) 

and that leads to 

6[ = -U(t) + ^s)sm\br + 6i) (C.36) 


<5^ = -7(0-l'3)sin2(6r + (52) (C.37) 


Returning to Eg. ( p .221 ) we find it reduces to 

= -acot(6r + 5)tanL). (C.38) 
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and combining that with Eq. (|C.32| ) yields 



[- sin(6r + 6) cos(D)<li + 26' csc(6r + 6) sin D] 
- cot{br + 6) tanD[^ cos(6r + 6) sin(D)l>i] 



so that 



e' = ^^i{tanDcos^{br + 6) - sin^(6r + 6) cot{D)) 



1 



-$i(sin^ D cos^(6r + 6) - sin^(6r + S) cos^(D)) 



2 sin D cos (D) 6 

^ l>i sin(6r + 6 + D) sin(6r + 6 - D) 



6 sin 2D 

From the definintion of $ we see that 



(C.39) 



cose sine \ / $3 $i \ / cose —sine 
— sine cose / \ $i —$3 / \ sine cose 

$3 cos2e + $1 sin2e -$3 sin2e + $1 cos2e 
-$3 sin 2e + $1 cos 2e -$3 cos 2e - $1 sin 2e 





' $3 ^'l ] 




I $1 -1.3 ) 



(C.40) 



So from this and Eqs. (|C.36| ) and Eg. ( |C.37|) we obtain the phase shift equations 
Eqs. ([5:251 , ^:2^ ) given in the text while Eq.( |C:39| ) gives us Eq.QO^). 



APPENDIX D: PHASE SHIFT EQUATION WITH THE COULOMB POTEN- 



TIAL 



We review here the necessary modification of our phase equations when we consider 
pp scattering [|16], When we study pp scattering, we must consider the influence of 
the Coulomb potential. The general form of the uncoupled Schrodinger-like equation with 
Coulomb potential is ||22| 



/(/ + 1) 2e^a 



+ A$ 



(D.l) 



where A<l> consists of the short range parts of the effective potential, a is the fine struc- 
ture constant. (Compare the Coulomb term with the first term on the right hand sides of 
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Eqs.( |2.27| , p.29| )). Due to the long range behavior of the potential in above equation, the 
asymptotic behavior of the wave function is 



u{r) const ■ sin(6r — ?7ln26r + A), (D.2) 

in which 

A = 5, + a;-y, (D.3) 

where ai = argr(Z + 1 + irj) is the Coulomb phase shift, here rj = — 

We describe here the variable phase method to calculate the phase shift with the Coulomb 
potential. Consider the two differential equations 

u" + {b^-W- W)u = 0, (D.4) 

and 

u" + {b"^ -W)u = 0, i = l,2 (D.5) 
in which u{0) = ui{0) = 0. Let 



— ^ _2e^^ 
r 



so that 

Mi(r) '^^^ const ■ sin(6r — ?7ln26r + A), 



U2i'i") '-^^ const ■ cos(6r — ?7ln26r + A), (D.7) 

where A = (Tq. 

Just as in the variable phase method, we obtain a nonlinear first order differential equation 

for the phase shift function 6i{r) such that Si{oc) = 6i, and 5/(0) = 0. This is done by 
rewriting u{r) as 

■u(r) = Q;(r) [cos 7 (r) Ml (r) + sin7(r)M2('")] (D-8) 

so that 

A = A + 7(cx)). (D.9) 
83 



Since we have rewritten u{r) in two arbitrary functions, we are free to impose a condition 
on u{r) 



u\r) = a' (r) [cos J {r)u[{r) + sin7(r)M2(r)]. (D.IO) 
Combining u{r) and u'{r) leads to 

where 7(0) = 0, and ui{r) — Fo{r],br) and U2{r) — Go{r],br) are the two Coulomb wave 
functions. With the Wronskian FqGq— F^Gq — h, we obtain, by differentiating, the differen- 
tial equation 

7'(r) = -iy(r)[cos7(r)Fo(?7,6r) + sin7(r)G'o(?7, 6r)]V&- (D.12) 

Note that for 

^ ^ A(A + 1) A(A + 1) ^ + _ 

^2 ' ^2 ^2 ^2 



Go(r?,6r)^-^, (D.13) 



we obtain the relation 



Co' 



7'(0) = -T^. (D.14) 



Letting 



A(A + 1) 



7(r) =/3(r) +r/(r), (D.15) 



where (5{r) is defined as 

(3'{r) = -^^^4^[cos7(r)Fo(77, 6r) + sin7(r)G'o(?7, br)f /h (D.16) 
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P{r) has the exact solution 



7(r) 



tan ^[ 



Fi{r], br)F^{r], br) - Fl{r], br)Fo{r], hr) 
Fi{r], br)G'f^{r], br) - Fl{r], br)Go{r], br) 



(D.17) 



with /3(0) = and p'{0) 



1(1+1} 



and P{oo) = a; — y — o"o, lead to 



6i = r]{oo). 



(D.18) 



Thus , if we solve 



+ A<l>] [cos(/3(r) + r]{r))Fo{r], br) + sin(/5(r) + r]{r))Go{r], br)f/b 



+ 




[cos/9(r)Fo(r7, br) + sin P{r)Go(r], br)f' /b 



(D.19) 



with the condition 77(0) = 0, we obtain the additional phase shift(above the Coulomb phase 
shift ) by integration to 77(00). 

There is no Coulomb scattering for the coupled triplet ^5*1 and ^Di states as a consider- 
ation of Pauli principal would show. 
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